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We consider orbital-tidal coupling in a binary stellar system to post-l-Newtonian order. We derive 
the orbital equations of motion for bodies with spins and mass quadrupole moments and show that 
they conserve the total linear momentum of the binary. Momentum conservation also allows us 
to specialize our analysis to the system's center-of-mass-energy frame; we find the binary's relative 
equation of motion in this frame and also present a generalized Lagrangian from which it can be 
derived. We then specialize to the case in which the quadrupole moment is adiabatically induced by 
the tidal field. We show how the adiabatic dynamics for the quadrupole can be incorporated into 
our action principle and present the simplified orbital equations of motion and conserved energy for 
the adiabatic case. These results are relevant to the gravitational wave signal of inspiralling binary 
neutron stars. 
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I. INTRODUCTION AND SUMMARY 
A. Background and motivation 

Inspiralling and coalescing compact binaries present 
one of the most promising sources for ground-based grav- 
itational wave (GW) detectors [l|. A primary goal in the 
measurement of GW signals from neutron star-neutron 
star (NSNS) and black hole-neutron star (BHNS) bina- 
ries is to probe the neutron star matter equation of state 
(EoS), which is currently only loosely constrained by elec- 
tromagnetic observations in the relevant denisty range 
p ~ 2-8 xlO 14 g cm 3 [2]. The EoS will leave its imprint 
on the GW signal via the effects of tidal coupling, as the 
neutron star is distorted by the non-uniform field of its 
companion. Many recent studies into the effects of the 
neutron star EoS on binary GW signals have been based 
on numerical simulations of the fully relativistic hydro- 
dynamical evolution of NSNS and BHNS binaries (see 
e.g. the reviews These simulations have largely fo- 

cused on the binaries' last few orbits and merger, at GW 
frequencies >500 Hz, and have investigated constraining 
neutron star structure via (for example) the GW energy 
spectrum @, effective cutoff frequencies at merger 
and tidal disruption signals [7] . 

As recently investigated by Flanagan and Hinderer Q , 
neutron star internal structure may also have a measur- 
able influence on the GW signal from the earlier inspiral 
stage of a binary's orbital evolution, at GW frequencies 
<500 Hz. While tidal coupling will produce only a small 
perturbation to the GW signal in this low-frequency adi- 
abatic regime, the tidal signal should be relatively clean, 
depending (at leading order) on a single parameter per- 
taining to the neutron star structure. This tidal de- 
formability parameter A is the proportionality constant 
between the applied tidal field and the star's induced 
quadrupole moment and is sensitive to the neutron star 
EoS. The measurement scheme proposed in Rcf. |8j is 
based on an analytical model for the tidal contribution 
to the GW signal, giving a linear perturbation to the 
GW phase proportional to A. At GW frequencies <400 
Hz, the model should be sufficiently accurate to constrain 



A to ^10%, with the largest source of error being post- 
l-Newtonian (1PN) corrections to the tidal-orbital cou- 
pling 8, 9]. It is our goal in this paper to lay the foun- 
dations for computing the tidal signal to 1PN accuracy. 

The modeling of tidal effects in GW signals from neu- 
tron star binaries can be divided into three separate prob- 
lems: (i) to calculate the deformation response of each 
neutron star to the tidal field generated by its companion, 
(ii) to calculate the influence of the tidal deformations on 
the system's orbital dynamics, and (iii) to calculate the 
gravitational waveform emitted by the system and incor- 
porate the corresponding radiation reaction effects in the 
orbital dynamics. This paper is focused on solving prob- 
lem (ii) to 1PN accuracy, while we leave problem (iii) to 
future work. A solution to problem (i) that is sufficient 
for our purposes has been given in Refs. [^-[l2j]. 

Tidal Deformability Computations: While computing 
tidal deformations of stars is a well studied problem in 
Newtonian gravity (l3| , it has recently been re-examined 
in the context of fully relativistic stars by Hinderer et 
al. foL [To | , D amour and Nagar [ll|, and Binnington and 
Poisson 121 . These authors used fully relativistic mod- 
els for the star's interior, with various candidate equa- 
tions of state, to calculate the perturbations to the star's 
equilibrium configuration produced by a static external 
tidal field. In Refs. [t| [loj], the results were used to de- 
termine the (electric-type) quadrupolar tidal deformabil- 
ity A, while Refs. [ill, G3 extended the analysis and in- 
cluded all higher-multipolar electric-type response coeffi- 
cients (for the mass multipole moments) as well as their 
magnetic-type analogues (for the current multipole mo- 
ments). While these results concern a star's response to 
a static tidal field, they are still applicable for inspiralling 
binaries at sufficiently low orbital frequencies; a star will 
approximately maintain equilibrium with the instanta- 
neous tidal field if the field changes adiabatically. 

Orbital Dynamics including Quadrupole: In their 
treatment of tidal effects in binary GW signals in Ref . |8[ , 
Flangan and Hinderer used Newtonian gravity to treat 
(the conservative part of) the system's orbital dynamics; 
they estimated that 1PN corrections to the orbits would 
modify the calculated tidal signal by ~10%. A general 
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formalism for calculating 1PN corrections to the orbital 
dynamics of extended bodies has been developed in a se- 
ries of papers [H], EH by Damour, Soffel, and Xu (DSX) 
and later extended by Racine and Flanagan (ltij . In this 
paper, we apply that formalism to determine and analyze 
the explicit 1PN translational equations of motion for bi- 
nary systems with quadrupolar tidal interactions. While 
such equations of motion have previously been presented 
in Refs. [ljl Gil, we find that those results differ from 
ours. We also extend previous results by analyzing 1PN 
momentum conservation for the binary system (which 
serves as a strong consistency check for our equations 
of motion, and one not satisfied by the results [T3, [HI), 
by specializing the equations of motion to the system's 
center-of-mass-energy frame, and by formulating an ac- 
tion principle for the orbital dynamics. 

Gravitational Wave Emission: To calculate the emit- 
ted GW signal and radiation reaction effects, Flanagan 
and Hinderer [8| used the quadrupole radiation approx- 
imation and associated 2.5PN radiation reaction forces. 
To consistently generalize their calculation in Ref. @ to 
(relative) 1PN order, it will be necessary to compute the 
3.5PN corrections to the GW generation, in addition to 
the 1PN orbital corrections considered in this paper; we 
leave these calculations to future work. 

The issue of tidal effects in inspiral-stage NSNS binary 
GW signals has also been recently addressed by means 
of numerical relativity simulations in Refs. jl9l l2fj||. and 
analytically, within the context of the effective one body 
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(EOB) formalism, by Damour and Nagar [2l|. Ref. 
proposes a contribution to the EOB "radial potential" 
that should account for the same 1PN corrections to the 
quadrupole-tidal interaction considered in this paper. In 
Appendix \K\ we demonstrate the equivalence of the EOB 
Hamiltonian of Ref. [2l| and our Lagrangian (|1.12j) for 
circular orbits; we also show how the EOB Hamiltonian 
can be extended to include 1PN tidal effects for non- 
circular orbits. 



B. Newtonian tidal coupling 

Gravitational tidal coupling arises from the interac- 
tion of the non-spherical components of a body's matter 
distribution with a non-uniform gravitational field. The 
non-sphericity is characterized at leading order by the 
body's mass quadrupole moment, 
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with p{t, x) being the mass density and x l {t) = x l — z J (i) 
being the displacement from the body's center of mass 
position x % = z*(t). The quadrupole is (at most) on the 
order of Q 1 ^ ~ MR 2 , with M being the body's mass and 
R its radius. Higher-order deformations are described 
by the octupole, Q ljk ~ J d 3 x px l x 3 x k ~ MR 3 , and 
higher-order multipole moments. The non-uniform field 



can be characterized by derivatives of an external Newto- 
nian potential 4> ex t(t, x); in a binary system, the leading- 
order potential is <fi cx t — —GM/r, where G is Newton's 
constant, M is the mass of the companion, and r the 
distance between the body and its companion. 

The non-uniform field of the companion produces tidal 
forces on the non-spherical body (in addition to the usual 
1 /r 2 force) according to 
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where M is the mass of the body or its companion, as- 
sumed here to be roughly equal, and the derivatives of 
the external potential are evaluated at the body's center 
of mass, x l = z % it). The contributions to the net force 
from the quadrupole and higher-order multipoles are thus 
seen to take the form of an expansion in R/r, the ratio 
of the size of the body to the orbital separation. When 
this finite-size parameter is small, as in the early stages 
of binary inspiral, the tidal force is well approximated 
by the quadrupolar term alone. A more detailed account 
of Newtonian tidal forces and related results is given in 
Sec. |H] below. 

In neutron star binaries, a quadrupole is induced by 
differential forces resulting from the non-uniform field of 
the companion. In the adiabatic limit, when the response 
time scale of the body is much less than the time scale 
on which the tidal field changes, the induced quadrupole 
will be given (to linear order in the tidal field) by 



'(<) = -\didj(t> ext (t,z), 



(1.3) 



with the constant A being the tidal deformability. This is 
related to the more often used dimensionless Love num- 
ber &2 by A = 2fc2-R 5 /3G, where R is the body's radius 

m 

Using Newtonian gravity to describe the orbital dy- 
namics and the adiabatic approximation to model the 
stars' induced quadrupoles, Flanagan and Hinderer Q 
calculated the effect of tidal interactions on the phase 
of the gravitational waveform emitted by an inspiralling 
neutron star binary and analyzed the measurability of the 
tidal effects. They found that Advanced LIGO should be 
able to constrain the neutron stars' tidal deformability 
to A < (2.0 x 10 37 gcm 2 s 2 )(D/50Mpc) with 90% con- 
fidence, for a binary of two 1.4 Mq neutron stars at a 
distance D from the detector, using only the portion of 
the signal with GW frequencies less than 400 Hz. The 
calculations of A for a 1.4 Mq neutron star in Refs. [9T4l2l]. 
using several different equations of state, give values in 
the range 0.03-1.0xl0 37 gcm 2 s 2 , so that nearby events 
may allow Advanced LIGO to place useful constraints on 
candidate equations of state. 
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Refs. 0, Q estimate the fractional corrections to 
the tidal signal due to several effects neglected by the 
model of the GW phasing used in Ref. [H, namely, 
non-adiabaticity (<1%), higher-multipolar tidal cou- 
pling (<0.7%), nonlinear hydrodynamic effects (<0.1%), 
spin effects (<0.3%), nonlinear response to the tidal 
field (<3%), viscous dissipation (negligible), and post- 
Newtonian effects (<10%). The largest corrections, from 
post-Newtonian effects in the orbital dynamics and GW 
emission, will depend on the neutron star physics only 
through the same tidal deformability parameter A used 
in the Newtonian treatment and thus can be easily in- 
corporated into the data anaysis methods outlined in 
Refs. EM- 



C. Post-l-Newtonian corrections 

For insprialling neutron star binaries with a total mass 
of ^3M Q at orbital frequencies of ^200 Hz (GW frequen- 
cies of ^400 Hz), the post-Newtonian expansion param- 
eter v 2 /c 2 ~ GM/c 2 r is ~0.1, so that the 1PN approxi- 
mation is well suited to describing relativistic corrections 
to the binary orbit. As discussed in depth in Sec. IHII bc- 
low, the 1PN orbital dynamics of a binary system with 
tidal coupling can be described by translational equa- 
tions of motion (EoMs) similar in form to the Newto- 
nian equations schematically represented in (|1.2p . giving 
the center-of-mass acceleration of each constituent body 
in terms of their positions and multipole moments. The 
1PN equations of motion add order 1 /c 2 correction terms 
to (|1.2[) which depend not only on the bodies' positions 
and mass multipole moments, M, Q Ij7c , etc., but 
also on their velocities and current multipole moments, 
S l , , etc. Expanding in both the post-Newtonian pa- 
rameter v 2 / c 2 and the finite size parameter i?/r, the 1PN 
equations of motion can be written schematically as 
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In the top line, the first term gives the usual point- 
particle (monopole) force of Newtonian gravity, and the 
second term represents its 1PN corrections. The last 
term of the top line denotes 2PN and higher post- 
Newtonian order corrections, which we neglect in this 
paper. Point-particle EoMs are in fact currently known 
up through post-3-Newtonian order [22j . 

In the second line of Eq. (|1.4[) , we have first the New- 
tonian quadrupole-tidal term, followed by its 1PN cor- 
rections (which are the subject of this paper), and finally 
contributions from the octupole and higher mass multi- 
poles (and their post-Newtonian corrections) which are 



surpressed by higher powers of R/r and which we neglect 
in our analysis [lq ]. 

The first term in the bottom line represents the 1PN 
spin-orbit coupling [15]. Though it would be convenient 
to be able to specialize to the case of non-spinning bod- 
ies when studying tidal interactions, tidal fields will in 
general lead to tidal torques on non-spherical bodies; a 
consistent 1PN treatment of tidal coupling thus requires 
the inclusion of 1PN spin-orbit coupling terms. The fi- 
nal term of the bottom line denotes 1PN contributions 
to the EoMs from the bodies' current quadrupoles and 
higher current multipoles, which we will not include in 
our analysis. 

The DSX treatment of 1PN celestial mechanics [i~l li~5l] 
provides a framework for calculating the orbital EoMs for 
bodies with arbitrarily high-order mass and current mul- 
tipole moments. In Ref. [15| . DSX applied their formal- 
ism to rederive the explicit 1PN EoMs for bodies with 
mass monopoles and current dipoles (spins). The calcu- 
lation was then extended to include the effects of the bod- 
ies' mass quadrupoles by Xu et al [13, EH- Racine and 
Flanagan (RF) [l6| later reworked the DSX formalism 
and presented explicit 1PN EoMs for bodies with arbi- 
trarily high-order multipoles, which can be specialized to 
the case of bodies with only spins and mass quadrupoles. 
We have found, however, that the final results of Xu et al 
and those of RF are in disagreement with each other and 
with our recent calculations. Some typos and omissions 
leading to errors in RF have been identified and are out- 
lined in an upcoming erratum; the results given in this 
paper agree with the corrected results of RF. In Sec. IIIII 
below, we review the essential ideas of the DSX formal- 
ism, following the notations and conventions of RF, and 
we outline the full procedure by which our results for the 
1PN EoMs are derived. 



D. Summary of results 

1. The M1-M2-S2-Q2 system 

Our results concern the 1PN gravitational interactions 
in a system of two bodies, which we label "1" and "2" . 
We model body 1 as an effective point particle, with 
a mass monopole moment only, while we take body 2 
to have additionally a spin and a mass quadrupole mo- 
ment. We consistently work to linear order in the spin 
and quadrupole, and our results can thus be easily gener- 
alized to the case of two spinning, deformable bodies by 
interchanging particle labels. We initially assume noth- 
ing about the bodies' internal structure or dynamics. Our 
primary assumption is the validity of the 1PN approxima- 
tion to general relativity in a vacuum region surrounding 
the bodies. 

The system's orbital dynamics can be formulated 
in terms of the bodies' center-of-mass worldlines z\(t) 
and z l 2 {t) and their multipole moments: the mass 
monopoles Mi(t) and M2(t), the spin S^t), an< ^ the mass 
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quadrupole Q 2 (t). 

The wordlines x l — z\(t) and x % = z 2 (t) parametrize 
the bodies' positions to 1PN accruacy in a (confor- 
mally Cartesian and harmonic) 'global' coordinate sys- 
tem (t, x l ), which tends to an inertial coordinate system 
in Minkowski spacetime as |a;| — > oo. The global coor- 
dinates and center-of-mass worldlines are defined more 
precisely in Sec. IIII Dl We use the following notation for 
the relative position and velocity: 



4> r=\z\ = y/SiJz^, 



z l /r, 



6 2 , 



J ll 



with dots denoting derivatives with respect to the global 
frame time coordinate t. 

The multipole moments— M x (f) , M 2 (t), S^(t), and 
Q 2 (i) for our truncated system — are defined in 
Sees. IIII Bl and HTTEl via a multipole expansion of the 1PN 
metric in a vacuum region surrounding each body fl(| . In 
the case of weakly self-gravitating bodies, these moments 
can be defined as intergrals of the stress-energy tensor 
over the volume of the bodies, as in Eqs. (|3 . 1 2[) : these def- 
initions coincide with those of the Blanchet-Damour mul- 
tipole moments introduced in Ref. [23| and used by DSX 

[UGl. The 

mass multipole moments (like Mi, M 2 , and 
Q 2 J ) are defined with 1PN accuracy, while the current 
multipole moments (like S 2 ) appear only in lPN-order 
terms and thus need only be defined with Newtonian ac- 
curacy. We will often denote the spin and quadrupole of 
body 2 by Si and Q y , dropping the "2" labels. 



2. General equations of motion and orbital Lagrangian 

The equations of motion for the monopoles Mi (t) and 
M2(t), the spin S l (t), and the worldlines z\(t) and z 2 (t) 
are determined by Einstein's equations alone, while that 
for the quadrupole Q % i (t) will depend on the details of 
body 2's internal dynamics and can initially be left un- 
specified. The mass monopole of body 1, the effective 
point particle, is found to be conserved to 1PN order: 

Mi = 0(c- 4 ), (1.5) 

while that of body 2 is not. As discussed in Sec. IIII Fl 
the lPN-accurate mass monopole M 2 can be decomposed 
according to 

M 2 = "M 2 + cr 2 [Ef + 3U Q ) + 0{c- 4 ). (1.6a) 

Here, "M 2 is the Newtonian-order (rest mass) contribu- 
tion, which is conserved: 

n M 2 = 0. (1.6b) 

The 1PN contributions to (ll.6a|) involve the Newtonian 
potential energy of the quadrupole-tidal interaction, 

3Mi 



2r 3 



(1.6c) 



and the Newtonian internal energy of body 2, E 2 nt , whose 
evolution is governed by the rate at which the tidal field 
does work on the body [24j |: 

3M A 



2r 3 
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(1.6d) 



(c.f. Sec. Ill F|) . Eqs. (jl.6[> ensure that M 2 satisfies the 
1PN evolution equation p.26!3.27p . The decomposition 
of the monopole M 2 in (|1.6a|) is essential for properly 
formulating an action principle for the orbital dynam- 
ics. The evolution of the spin S % is determined by the 
(Newtonian-order) tidal torque formula: 



S l = 



3Mi 



% i k Q3 a n a n k + 0{c~ 2 ), (1.7) 



as in (|2.45|) . Finally, the 1PN translational equations of 
motion, which govern the evolution of the worldlines z\ 



and z\, are of the form 



z\ = T[(zi ,v{,vi,Mi,M 2 ,S> ,Q> k ,& k ,& k ), 
z\ = r 2 (z^v{,4,M 1 ,M 2 ,S j ,Q jk ,Q :ik ), 



and are given explicitly by Eqs. (|3.3ip . 

In Sec. IIV1 we define and calculate the lPN-accurate 
mass dipole moment of the entire system M* ys (t), given 



in Eq. 



We find that the condition M, 



sys 



0(0, 



required by Einstein's equations and reflecting the con- 
servation of the system's total momentum, is satisfied 
as a consequence of the orbital EoMs (|3.3ip : this serves 
as a non-trivial consistency check for our results. The 
conservation of momentum also allows us to specialize 
the EoMs to the system's center-of-mass(-energy) (CoM) 
frame, which can be defined by the condition Mj! ys = 
as in Sec. IVAl The two EoMs (|3T3T1) for the worldlines 
z\ and z 2 can then be traded for the single EoM for 
the CoM- frame relative position z % — z\— z{, given in 
Eq. dHJ. 

Our results can be most compactly summarized by giv- 
ing a Lagrangian formulation of the CoM-frame orbital 
dynamics, as discussed in Sec. (TV Bp . We find that the 
CoM-frame orbital EoM (|5 .0[) can be derived from the 
Euler-Lagrange equation 



d _ d d 
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with a generalized Lagrangian £ or b given by 

£orb = + £-S + £q, 

with the monopole part, 



(1.8) 



(1.9a) 
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+ 0(c- 4 ), 



the spin part, 
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+ (9(c- 4 ), (1.9c) 



and the quadrupole part, 
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We use here the notation 

M = M x + n M 2 , xi = Mx/M, X2 = n M 2 /M, 
fx = Mx n M 2 /M, 77 = X iX2 = n/M, 

with M being the total (conserved) Newtonian rest mass, 
/i the reduced mass, and r\ the symmetric mass ratio. The 
dimensionless coefficients A±-Ag appearing in (jl.9dj) are 
functions only of the mass ratios xi an d Xi given by 

oa . 

3. Adiabatic approximation for the induced quadrupole 

The above results concerning the orbital EoM and its 
Lagrangian formulation are valid regardless of the inter- 
nal structure of body 2, i.e. for arbitrary evolution of its 
quadrupole Q l ^{t). In Sec. IVI[ we discuss a simple adi- 
abatic model for the evolution of Q lJ . In the adiabatic 
limit, the quadrupole responds to the instantaneous tidal 
field according to 

Q i \t) = \G 2 j {t), (1.10) 

where G 2 3 {t): given by Eq. (|6.6bj) . is the quadrupolar 
gravito-electric tidal moment of body 2, a 1PN general- 



ization of the derivatives of the Newtonian potential in 
Eq. (|1.3p . and A is the lPN-accurate tidal deformability. 
With the quadrupole given by (|1.10[) . the tidal torque 
(fl~7]) vanishes (c.f. Eq. (f!T46]> ). so that the spin S i is 
constant; thus, in the adiabatic limit, we can specialize 
to the case S l = without generating inconsistencies. In 
Sec. IVI A) we show that the adiabatic evolution for the 
quadrupole (|1.10j) can be derived from a Lagrangian that 
adds to the orbital Lagrangian an internal elastic poten- 
tial energy term which is quadratic in the lPN-accurate 
quadrupole: 



C = C mh [z\ Q«] - ^-Q ab Q ab + 0(c- 4 ), (1.11) 



with £ or b given by Eq. (|1.9p with S l = 0, and with 
E 2 nt = (l/4X)Q ab Q ab + 0(c~ 2 ). (Note that any addi- 
tional constant contribution to the internal energy E 2 nt 
can be included as a 1PN contribution to the constant 
n M 2 in Eq. H~6al ) Substituting the solution ([i~T0]) for 
Q lJ (t) into this Lagrangian, we obtain a simplified La- 
grangian involving only the CoM-frame relative position 
z*{t): 



,M^ + ^(i + A) + ii 



hv H 

r 



, 2+& A) + M(„ a+& A 



(1.12) 



r 



with A = (3%i/2%2)A, and with the dimensionless coeffi- 
cients 

9 = (l-3ry)/8, 

6x = (3 + 77)/2, 

02 = fi/2, 

03 = -1/2, 

6 = (xi/2)(5 + X 2 ), 

6 = -3(l-6 X 2 + xi), 

6 = -7 + 5x2- (1.13) 



The orbital EoM resulting from this Lagranigan, which 
can also be found by substituting the adiabatic solu- 
tion (jl.lOp for Q* J '(t) directly into the general orbital 
EoM ([53]) . is given by Eq. ([STCP]) . The conserved en- 
ergy E(z, v) derived from the Lagrangian (|1.12l) is given 
by Eq. (I6.13[) . In the case of circular orbits, we find the 
gauge-invariant energy-frequency relationship 
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E{u) = p(Afw) 2/3 
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This result, along with others from this paper, will be 
useful in calculating the phasing of GW signals from in- 
spiralling neutron star binaries. 



E. Notation and Conventions 

We use units where Newton's constant is G = 1, but 
retain factors of the speed of light c, with 1/c 2 serv- 
ing as the formal expansion parameter for the post- 
Newtonian expansion. We use lowercase latin letters 
a, b, . . . for indices of (three-dimensional) spatial ten- 
sors. Spatial indices are contracted with the Eucliden 
metric, v l uf = 8ijV*w J , with up or down placement of 
the indices having no meaning. We use uppercase latin 
letters to denote multi-indices: L denotes the I indices 
a\ai ... a/, K denotes the k indices 6162 • etc. For a 
given vector v* or for the partial derivative operator di, 
we use multi-indices or explicit sequences of indices to 
denote their tensorial powers: 



V L = v a ia2 ... ai =v a lv a 2 __ y a K 



d K = d, 



■ 0, 



b, ■ 



(1.15) 



and, for example, u y = i/?;- 7 . We also use v 2 = v n 
and V 2 = da. Multi-indices are also used for sets of 
distinct tensors of varying rank, {M, M a , M ab , . . .}, with 
M L = M aiCL2 '" a > denoting the tensor of rank I. We use 
the Einstein summation convention for both individual 
indices and multi-indices. Derivatives with respect to a 
time coordinate t are denoted by dt or by overdots. 

We use angular brackets to denote the symmetric, 
trace- free (STF) projection of tensors 14j: 



rj-i<ab> 
rp<abc> 
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(1.16) 



and so on, with parentheses denoting the symmetric pro- 
jection. For a STF tensor S L — S <L> and general tensor 
T L , note that S L T L = S L T <L> . 



II. NEWTONIAN TIDAL INTERACTIONS 



bodies. We consider in particular the case of a binary 
system containing a point particle (body 1) and an ex- 
tended deformable star (body 2), working to quadrupolar 
order in the star's multipole series. We also discuss an 
action principle formulation of the orbital dynamics, the 
process of energy transfer between the gravitational field 
and the deformable body, and the evolution of the body's 
spin due to tidal torques. Finally, we discuss the evolu- 
tion of the body's tidally induced quadrupole moment in 
the adiabatic limit. 



A. Field equations 

In Newtonian physics, the scalar potential <p(t, x) 
obeys the Poisson equation, 



(2.1) 



with p(t, x) being the rest mass density of matter. The 
influence of the gravitational field on matter is described 
by the test particle acceleration x l = —di<p, or more gen- 
erally, by Euler's equation supplimented by the continu- 
ity equation (the conservation of mass), 

d t (pv l ) + djipvV + V) = -pdi4>, (2.2) 



p + d i (pv i )=0, (2.3) 

with v l (t,x) being the matter's velocity field, and 
i IJ (i, x) the material stress tensor. Together, Eqs. (|2.U 
12. 3|) provide a complete description of Newtonian grav- 
itational interactions. However, they do not in general 
form a closed set of evolution equations for the fields (j>, 
p, v % , and £ y ; one needs also to specify the matter's inter- 
nal dynamics, in particular concerning the stress tensor 
tij. (In the simplest cases, one can fix t u by an algebraic 
equation, like i iJ = for 'dust' or i IJ = pS^ with p(p) 
being the pressure for an isentropic perfect fluid.) Still, 
one can derive many useful results, like the form of the 
translational equations of motion for a system of gravitat- 
ing bodies, while leaving the matter's internal dynamics 
unspecified. 



In this section, we review the standard treatment of 
tidal coupling in Newtonian theory [l3| : the post-1- 
Newtonian treatment given subsequent sections makes 
extensive use of these Newtonian-order results. We define 
the multipole moments and tidal moments of an extended 
object and use them to derive the orbital (or transla- 
tional) equations of motion for systems of gravitating 



B. Multipole moments 

We consider iV isolated celestial bodies, i.e. regions of 
space containing matter (p ^ 0) surrounded by regions 
of vacuum (p — 0), and label the bodies by an index A, 
with 1 < A < N. The potential that is locally generated 
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by body A, which we will call the internal potential A ' ul 



A i 



is given by the standard solution to (|2.ip as an integral 
over the volume of the body: 



j3m/ p(t,x') 



IX — X' 



(2.4) 



We can express this potential as a multipole series around 
a (moving) point x l = z A (t) by using the Taylor series 



i 



\x — X 



1=0 
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\x - Z A \ 



(2.5) 



where L = a\a2 ... a; is a spatial multi-index, denoting 
here tensorial powers of the vector (x' — z A ) z and of the 
operator di (c.f. Eq. (|1.15|l ). Using the Taylor series 
(12.5[) in (|2.4p allows us to write the internal potential, for 
points x 1 exterior to the body, in the form 



fc (^) = -E 



1=0 



\x - z A (t)\ ' 



(2.6) 



about which the multipole expansion is centered to be 
fixed to the body's center of mass position: 



1 

~M A 



d 3 x p{t, x) x l 



(2.10) 



which is equivalent to (|2.9b|) . Next comes the I = 2 
term involving the body's mass quadrupole tensor Q l A (t) 
(|2.9c|) . which we have renamed M A —> Q A to accord 
with common convention. The higher-order terms in the 
multipole series are supressed by increasing powers of 
the finite size parameter R/\x\, where R is the size of the 
body. 

We note that the continuity equation (|2.3p implies the 
constancy of the body's total mass, M A = 0. The Eulcr 
equation (|2 . 2[) similarly constrains the evolution of the 
mass dipole MjL thus determining the body's transla- 
tional equation of motion, as discussed in Sec. Ill Dl below. 
Eqs. (I2.2[) and (I2.3|) do not, however, fully determine the 
evolution of the quadrupole and higher multipoles, which 
will depend on the body's internal dynamics. 



with 



Tidal moments 



M%(t) = / d 3 x p(t, x) [x ~ z A (t) 



<L> 



(2.7) 



The quantities M A are the mass multipole moments of 
the body about the worldline z A (t). They are symmet- 
ric, trace- free spatial tensors of varying order I, M A — 
M a ia2 ...a,_ The STp property f n ows f rom t he fact that 

dj,|a;| -1 is an STF tensor, because partial derivatives 



commute, and because c?,-,-|a;| 



V 2 | 



0. As 



the multipole moments M A are contracted with the STF 
tensors &l\x — z A \~ 1 in (|2.6p . only their STF parts will 
contribute to the potential (f> A t ; this is why the STF pro- 
jection (denoted by angular brackets) has been included 
in the definition of the multipole moments in (I2.7p . 

The leading-order terms in the multipole series (12. 6P 
can be written more explicitly as 



Ma 
' \x\ 



Qi_dij 



1 



+ 



(2.8) 



Having described the potential generated by an iso- 
lated body with its multipole moments, we can similarly 
describe the potential felt by the body with tidal mo- 
ments. Given a collection of several bodies indexed by 
B, each giving rise its own intrinsic potential of the form 
(I2.6p . we define the external potential felt by a given 
body A to be the sum of the potentials due to all the 
other bodies: 



B^A 



(2.11) 



The body's tidal moments 1 A (t) are then defined as 
coefficients in the Taylor expansion of the external po- 
tential about the center-of-mass position z\: 

oo 

x) = - J2 Tf Gt A {t) [x z A (t)] L , (2.12) 
i=o 



with 



M A 



= M\ = 



Qi=M i i = 



(2.9a) 
(2.9b) 



d*xp ( x l x J - ~\xy8^ ) , (2.9c) 



d x p, 



d x p x % 



and x l — x % — z\. First is the monopole term (I = 0), 
generated by the total mass of the body M A (|2.9a|) . and 
giving rise to a Coulomb- type potential in (12.81) . We have 



omitted the dipole term (I = 1) in (|2.8p because M\ can 
always be made to vanish by choosing the point z A (t) 



G^ A (t) = -d L cbT(t,x)\ 



5a(*) ' 



(2.13) 



1 The subscript g, standing for 'global,' has been included here 
to avoid confusion with tidal moments introduced in our post- 
Newtonian treatment below. We introduce there a set of body- 
frame tidal moments G A , defined in an accelerated reference 
frame attached to the body, and a set of global-frame tidal mo- 
ments Gg A , defined in an (assymptotically) inertial frame. The 
tidal moments defined in I I2.13H coincide with the latter at New- 
tonian order. While we have chosen to work exclusively in an 
inertial frame in our Newtonian treatment here, an analogous 
Newtonian treatment that uses accelerated frames can be found 
in Sec. Ill of Ref. [[!. 
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Like the multipole moments, the tidal moments are STF 
tensors, G g A = G< h A , as can be seen from the definition 

(|2.13|) and the fact that V 2 05f* = everywhere outside 
the bodies B ^ A. We see that G g:A is simply (minus) 
the potential at the body's center, and G z g A is the would- 
be test particle acceleration —dicj)^. For I > 2, the G gA 
are higher-order derivatives of the potential that will give 
rise to tidal forces on a nonspherical body. 

The tidal moments of a body A can be expressed in 
terms of the multipole moments of the other bodies B ^ 
A by combining 0Ji\) , (j2~TT1) . and ([2~T3|l . with the result 



EE' 



k\ 



"Mr [ , 



B U KL 



za — z B \ 



(2.14) 



find the following EoMs: 

M^l = M 1 M 2 d^—±—+ 1 -M 1 Qi k d^l ' 
\Z\ — z 2 \ l 



zi- z 2 



M 2 z\ = M 2 M l af ) • w.;i -' 1 



\z 2 - Zi\ 



z% - Zl\ 



Defining the radius and unit vector associated with the 
relative position, 

z' = 4-zi, r = \z\, n l = z l /r, (2.17) 

and using the general identity, 



1 n <L> 

^- = (-i)^-i)!<— 



these EoMs can be written as 



(2.18) 



where d\ A) = d/dz A , and d K L = d bl . . . d bk d ai ■ ■ ■ d a 



D. Translational equations of motion 

A primary advantage of the language of multipole and 
tidal moments is that it allows one to take the PDEs (12.11 
12. 3[) governing the evolution of the fields p, v l , , and 
4> and extract from them ODEs for the center-of-mass 
worldlines z l A (t) of a collection of gravitating bodies. To 
this end, we consider a body A with multipole moments 
M A defined by (|2.7|) , in the presence of an external poten- 
tial 05f* generated by other bodies B ^ A according to 
(j2~T2"|) and (|2~T4l . The body's translational EoM can be 
found by applying two time derivatives to the definition 
of z\ in (|2.10p , using the Euler and continuity equations 
(|2.3p and (|2.2p . and intergrating by parts. The result is 
an expression for the body's center-of-mass acceleration, 



M A z\ 



d 3 x p di 



(2.15) 



which can be rewritten in terms of the body's multipole 
and tidal moments by using (|2.12l) and (12.71) : 



M x z\ 



-M 2 i. 



2Z 2 

r 2 2r 4 



-71, 



(2.19) 



In this form, it is evident that the total momentum p l = 
Miz\+M 2 z 2 is conserved, and that these two EoMs for z\ 
and z\ can be traded for the EoM of the relative position 

z % = z\ -~ z\: 



2M 2 r 4 



(2.20) 



= - ^—Q ]k (5n^ k -2S^n k ) 

r 2 2M 2 r 4 ^ y ' 

where M = Mi + M 2 is the total mass. In the second 
line, we have used Eq. (|1.16[) to expand n <lJ and the 
fact that QJ k is STF. 

With Eq. (|2.20[) , we have reduced the description of 
the binary system's translational dynamics to a single 
ODE. Still, we can only solve this ODE if we know the 
time evolution of the quadrupole moment Q IJ (i), which 
requires a detailed model of body 2's interior. We will 
describe a simple adiabatic model for Q 1 ^ in Sec. Ill Gl 
below. 



M A z\ = 



T\ M aG\.a 

1=0 



(2.16) 



= M A Gl A + l -Q A k Gf A 



sA 



The first term in the second line represents the force that 
would act on a freely falling test mass M A , while the 
second term gives the leading-order tidal force. 

To render the EoM (|2.16j) fully explicit, one can use the 
expressions for the tidal moments (12. 14)) in terms of the 
multipole moments and worldlines of the other bodies. 
Considering the case of a two-body system A = 1,2, with 
body 1 having only a monopole moment M 1; and body 2 
having a monopole M 2 and a quadrupole Q 4 - 7 = Q 2 , we 



E. Action principle 

As for any Newtonian system, a Lagrangian for a col- 
lection of gravitating bodies can be constructed from 
C = T — U, with T being the kinetic energy and U the po- 
tential energy. The total kinetic energy receves separate 
contributions from each body, T — J2 A Ta, and each T A 
can be split into a contribution from the center-of-mass 
motion of the body and an internal contribution: 



T 



hit 



1 



1 



(fx p V 



1 



-M A z\ + T 



hit 

A i 



= o / d A xp{v- z A y 



(2.21) 
(2.22) 
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Here, we have used the fact that J A d 3 x p v l = Maz a , im- 
plied by the continuity equation (j2.3[) and the definition 
of the center-of-mass position z\ in (|2.10j) . The gravi- 
tational potential energy U = ^2 a ^ A can be similarly 
split into external and internal parts: 



d 3 x p<f> = uf t + uT, 



J rCXt 

U A 



d 3 x p 



1 

2 J A 

- I d'xp, 



(2.23) 
(2.24) 

oEn M i G g,A- (2-25) 



1=0 



In the last line, we have used (|2.12l) to express 



terms of the tidal moments and (I2.7|) for the definition of 
the multipole moments. 

While the system's total potential energy will also re- 
cieve non-gravitational contributions from the internal 
structure of each body, we can lump these contributions, 
along with T A nt and U A l as defined above, into an inter- 
nal Lagrangian C 1 ^ for each body. Wc can then write 
the total Lagrangian for an TV-body system as 



(We have omitted the internal Lagrangian for body 1 
as it is completely decoupled from the rest of the sys- 
tem.) Varying this action with respect to z\ and z\ 
leads to their EoMs found in (|2.19p . Alternately, varying 
with respect to their separation z l = z\ — z\ gives the 
relative EoM (|2.20[) . while the system's center of mass 
(M\z\ + M 2 z 2 )/M is found to be cyclic. Specializing 
the Lagrangian (|2.27[) to the center-of-mass frame, where 
M Y z\ + M 2 z l 2 = 0, gives 



C = 



flZ 



pM 



-Uq 



£int 
2 ' 



(2.29) 



with p = M\M 2 /M being the reduced mass. 

While leaving the functions C 2 nt (q 2 , 5f) an d Q lJ (12) 
unspecified, we can still use the Lagrangian (|2.27[) to 
write down EoMs for body 2's internal configuration vari- 
ables q 2 : 



d dC 2 nt 



dC 2 nt 



dt dqg dq% 



-G l - 



dQ l - 



(2.30) 



which will be useful in the next subsection. 



c = E ( \ M ^A 



1 



1 



2^l\ 

1=0 



Of 



(2.26) 



The bodies' center-of-mass worldlines z l A {t) enter this La- 
grangian through the translational kinetic energy terms 
and through the external gravitational potential energy 
terms (via the tidal moments); the internal Lagrangians 
C 1 ^, however, are independent of the worldlines z\, by 
construction. The will be functions of some set of 
internal configuration variables q A (and their time deriva- 
tives) for each body, which will include e.g. Euler angles 
for the orientation of the body, vibrational mode ampli- 
tudes, etc., depending on the model of the body's internal 
structure. (In full generality, the proper internal configu- 
ration variables are the fields p, v % and i lJ , subject to 
(|2.10p as a constraint.) The bodies' multipole moments 
M A , for I > 2, appearing in the gravitational potential 
energy terms in (|2.26p . will be functions of these same 
internal variables q^. Together, the z % A and q% for all 
bodies A, form a complete set of dynamical variables for 
the TV-body system. Varying the action S — J Cdt with 
respect to the worldlines z l A reproduces the translational 
EoMs (|2.16[) . Determining the evolution of the variables 
q A , and hence the moments M A for I > 2, will require a 
model for £ i £ t (q A [,q%) and M%(q%). 

Specializing to the two-body M\-M 2 -Q 2 case and using 
(prnjl . l[2"TT) . and ([2~T8]) . the Lagrangian ([2T26]) becomes 



1 9 1 9 MiM 2 



C/ Q +4 nt - (2-27) 



F. Energy 

Continuing to specialize to the two-body M\-M 2 -Q 2 
case, we can construct a conserved energy for the system 
from 

dC dC dC 

E = T +v=£% + W% + ^- £ > (2 - 31) 



with summation over a implied. Using the (CoM-frame) 
Lagrangian (|2.29l) . we find 



E 



pz 2 pM 



Uq + E? 



(2.32) 



where Uq is given by (|2.28[) . and the internal energy of 
body 2 is given by 



dC 2 nt 

dq% 



12 



£, nt 



(2.33) 



This internal energy will generally have several contribu- 
tions: internal gravitational potential energy, rotational 
kinetic energy, vibrational kinetic and potential energy, 
thermal energy, etc. Nonetheless, the rate at which en- 
ergy is exchanged between the interior of body 2 and its 
surroundings, via the gravitational tidal interaction, is 
a function only of the orbital separation, Mi , and Q 1 ^ . 
Differentiating (|2.33[) with respect to time and using the 
EoM (|2.30[) for the internal variables q 2 , we find 



where Uq is the potential energy of the quadrupole-tidal 
interaction: 



3Mi 



(2.28) 



E 2 at 



<12 



d dtif 



d£ 2 nt 



dt dq% dq% 
2 ^ ~ 2r 3 



■n 



dQ'< 



(2.34) 
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The energy transfer described by (I2.34[) is often referred 
to as tidal heating (see e.g. HI])- This expression for the 
power delivered to the body is valid (in the quadrupolar 
approximation) regardless of the body's internal dynam- 
ics. Using (|2~M| and the orbital EoM (pT20|) , one can 
confirm that the binary system's total energy (|2.32p is 
conserved. 



G. Adiabatic approximation 

While we have thus far left unspecified the internal 
dynamics for the deformable body 2, which determines 
the evolution of the quadrupole, we now specialize our 
analysis to the case where Q y (t) is adiabatically induced 
by the tidal field. This will lead to a closed system of 
evolution equations for the binary. In the adiabatic limit, 
when the body's internal dynamical time scales are much 
less than the orbital period, the qudrupole will respond 
to the instaneous tidal field according to 

Qv(t) = \G%(t), (2.35) 

where A is the tidal deformability, and Gg b 2 (i) is the tidal 
tensor given in Eq. (|2.28|) . As discussed in Sec. II Bl and in 
more detail in Refs. [8|, Q, the relation f|2 . 35[) should be 
valid to ~ 1% for neutron star binaries at GW frequencies 
<400 Hz. 

The adiabatic evolution of the quadrupole can be inco- 
prorated into our action principle (|2.29[) by taking Q 1 ^ (t) 
to be our lone internal configuration variable (g 2 ), and by 
taking the internal Lagrangian £ 2 nt to contain a simple 
quadratic potential energy cost for the quadrupole: 

C[z\ Q ij ] = + ^ - U Q + 4 nt (2-36) 

_ flZ 2 flM lfiab r^ab 1 nabnab 

- — + — + 2 G ^Q Q ■ 

Varying the action with respect to the orbital separation 
z l still gives the orbital EoM (|2.20p . and varying with 
respect to the quadrupole Q v reproduces the adiabatic 
evolution equation (|2.35[) . Using Eq. (|2.35[) to replace Q tj 
and Eq. (|2.28[) for Gg b 2 , the Lagrangian can can written 
solely in terms of z % as 

This Lagrangian leads to the orbital EoM (ET2"0")) with Q ij 
replaced by its adiabatic value (|2.35l) . 

which shows that the tidal coupling results in an attrac- 
tive force. For circular orbits, with a 1 — -~ruP"n l ', we find 
the radius-frequency relationship 

^^T^r^-ldfwj' (2 ' 39) 



to linear order in the tidal deformability A. 

The internal energy E 2 nt (|2.33l) . in the adiabatic ap- 
proximation, is given by 

E 'T = YrQ ab Q ab , (2-40) 
4A 

up to a constant, and satisfies the tidal heating equation 
([til]) by virtue of Eq. (|2~35|) . (We should note that there 
is actually no 'heating' taking place here, as this model 
neglects dissapative effects and is completely conserva- 
tive.) Using Eq. (|2.35l) . the binary system's total energy 
E (|2.32l) can be written as 

nz 2 uM ( 3AMi \ 

E = t-7T ~ — 1 + „ , 2.41 

2 r \ 2r 5 M 2 J 

in the adiabatic model, and is conserved by the orbital 
EoM (|2~38)) . Then, using i 2 = r 2 uj 2 and Eq. (12.39[) . we 
find 

for the circular-orbit energy-frequency relationship. 
H. Spin 

In anticipation of the 1PN treatment of the binary's 
orbital dynamics, in which a body's angular momentum 
(or spin) has a direct influence on the orbit, it will be 
useful to discuss the evolution of an extended body's spin 
at Newtonian order. The spin of a body A about its CoM 
worldline z\ is defined by 

S%(t) = e abc J d 3 xp(t,x)[x b - z b A (t)]v c (t lX ), (2.43) 

with p being the mass density and v c the velocity field. 
Taking a time derivative of this equation, using the Euler 
and continuity equations (|2.2p and (|2.3|) and the defini- 
tion of z\ in (|2.10[) . and integrating by parts, we find 

no abc I j3_ ( b b \ o jcxt 

GO 

= e abc J2-n M A LG £A- (2-44) 

1=0 

In the second line, we have used the definitions of M A and 
A in (12.71) and (|2.13l) . This formula gives the torque 
on the body due to tidal forces. As it is not directly 
relevant to our purposes, we will not discuss a Lagrangian 
formulation of the Newtonian rotational dynamics. 

Applying Eq. (|2.44l) to the M x -M 2 -Qi system, we find 
that the tidal torque on body 2 is given by 

ga = e abc Q bd G c^ (2 45) 

with the tidal tensor G c g d 2 given by (|2~2gj) . Eq. (EH51) is 
valid (in the quadrupolar approximation) regardless of 
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the internal dynamics of body 2. In the special case of 
an adiabatically induced quadrupole, as in Eq. (|2.35p , we 
find 



S" 



\ „abcr~ibd fHcd 



o, 



(2.46) 



so that the spin is conserved. 



III. POST-NEWTONIAN TIDAL 
INTERACTIONS 

A. Overview 

The Newtonian theory of gravity arises as a limit- 
ing case of general relativity (GR). In the limit of small 
source velocities and weak gravity, the spacetime metric 
of GR takes the form 



ds 2 



1 + ^ ) c 2 dt 2 + 5 %: >dx l dx ] + 0{c-' 2 ), (3.1) 



with <j)(t, x) being the Newtonian potential. This expres- 
sion represents a perturbation expansion of the theory 
with 1 /c 2 playing the role of a formal expansion param- 
eter. At leading order in 1/c 2 , Einstein's equation and 
covariant stress-energy conservation for the metric (I3.1[) 
reproduce the Poisson, Euler, and continuity equations 
(12. 1112. 3|) — the basic equations of Newtonian gravity. 

The first post-Newtonian (1PN) approximation to GR 
continues this perturbation expansion to next-to-leading 
order in 1/c 2 . The 1PN metric can be written as 



ds 2 = - 



20 



2C 

c 2 dt 2 + -\dtdx 1 



+ ( 1 - ~f J S^dx'dx 1 + 0(c~ 4 ), 



(3.2) 



(c.f. Weinberg 26]), with two new degrees of freedom 
appearing: a 1PN (three-)vector potential C{t, x) (often 
called the gravito-magentic potential), and a 1PN scalar 
potential ip(t, x). We will find it convenient to work with 
a single scalar potential $(t,x) which has and ip as 
its Newtonian- and lPN-order parts (and hence a hidden 
c-dependence) , 



4> 



+ c' 2 i> + 0(c- 4 ), 



(3.3) 



so that the metric can be written as 

ds 2 = - 



2$ 2$ 2 

1 + — + 



2C 

c 2 dt 2 + Ardtdx 1 
c z 



2$ N 

1 - — ] 8 l Hx l dx J + 0(c" 4 ). 



(3.4) 



We choose to work here in conformally Cartesian co- 
ordinates (see e.g. [3]), which is already implicit in the 
form (13.4[) of the metric, and to adopt the harmonic gauge 
condition: 

dn(V=ggn = o «*• 4* + diC = 0(c- 2 ). (3.5) 



In this gauge, one finds that Einstein's equation for the 
metric, at next-to-leading order in 1/c 2 , is equivalent to 
the following linear field equations for the potentials: 

V 2 $ = 4ttT u + c~ 2 (inT 11 + <t^ +0(c- 4 ), (3.6a) 

V 2 C = 167rT" + 0(c- 2 ), (3.6b) 

where T^ v are the contravariant components of the 
stress-energy tensor in the (t, x l ) coordinate system. 
Note that the T u component must include both 0(c°) 
Newtonian and 0(c~ 2 ) post- Newtonian contributions, 
while the components T tl and T y and the quantity $ 
are needed only to Newtonian order. The influence of 
the gravitational field on matter is governed by covari- 
ant stress-energy conservation: V^T^ = 0; the 1PN- 
expanded form of this equation can be found in Appendix 
DofRF[l6|. 

In the remainder of this section, we review the treat- 
ment of tidal interactions within this post-l-Newtonian 
framework. We employ a formalism, originally devel- 
oped by DSX [3, EH and later expounded upon by RF 
[16j . that uses multiple coordinate systems to describe 
the global motion and local structure of extended bod- 
ies. We attempt to present here the primary ingredients 
and broad logical flow of this formalism, which are essen- 
tial for properly interpreting the results stated in Sec. lIDI 
above. 

We begin in Sec. IIII Bl by presenting a general solution 
to the lPN-order Einstein equations (13.61) which gives 
the spacetime metric in a vacuum region surrounding an 
astronomical body A. The solution (I3.8[) is parametrized 
by and defines the body's multipole moments and tidal 
moments. The mass and current multipole moments M A 
and S A characterize the body's internal structure, and 
the gravito-electric and -magnetic tidal moments G A and 
H A characterize the external gravitational fields felt by 
the body. 

In Sec. IIII Cl we discuss the gauge freedom in the 1PN 
metric, summarized by the parametrization of a general 
1PN coordinate transformation in Eq. p. 131) . In 1PN 
celestial mechanics, it is advantageous to use and trans- 
form between two types of coordinate systems: global co- 
ordinates [t^x 1 ) used to describe the motion of multiple 
bodies, and body-adapted coordinates (sa,Ua) use cl in 
the local description of a given body A. We discuss how 
to fix all 1PN coordinate freedom in the body-adapted 
coordinates (sA,y\) by enforcing the body-frame gauge 
conditions (|3.15[) . The body- frame multipole moments 
M^(s A ) and S^[s ) — the moments defined by the mul- 
tipole expansions in Sec. IIII Bl using the body-adapted 
coordinates — then become unique and meaningful de- 
scriptors of a body's internal structure. 

In Sec. MIDI we discuss the form of the metric in the 
global coordinates (t,x l ). It is written in terms of a set 
of global- frame multipole moments M^ A (t) and Z^ L A (t) 
and tidal moments A (t) and Y^ L A (t) for each body A, 
which differ from the body-frame moments. The relation- 
ship between the global- and body-frame moments is de- 
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termined by the transformation f|3 . 13[) between the global 
and body-adapted coordinates; the moment transforma- 
tion formulae are presented in full detail in Appendix [5] 
The functions parametrizing the coordinate transforma- 
tion, or the worldline data T>a f|3 . 14[) , are seen to take on 
the role of configuration variables for body A with respect 
to the global frame. Among other things, they determine 
the body's center-of-mass worldline, x % = z\(t). 

In Sec. MI El we discuss the 1PN single-body laws of 
motion Q3.20P which govern the evolution of a body's 
mass monopole Ma<, mass dipole M A , and current dipole 
(or spin) S\. These laws reflect the conservation of en- 
ergy, momentum, and angular momentum and can be de- 
rived from stress-energy conservation at 1PN order [Tlj . 
or equivalently, from Einstein's equation at 2PN order 
16] . A body's translational equation of motion — an ODE 
for its global- frame center-of-mass worldline z A (t) — can 
be deduced from the law of motion for its body-frame 
mass dipole. The result is an expression for the accelera- 
tion z A , for each member A of an TV-body system, written 
in terms of the body-frame multipole moments M A and 
S A and global-frame worldlines z\ of all the bodies A. 

Finally, in Sec. IIIIF1 we specialize our discussion to 
the case of a two-body system with a body 1 having 
only a mass monopole Mi, and a body 2 having a mass 
monopole M%, a mass quadrupole Q 1 ^ = Q lJ , and a spin 
$2 = S l . We present and discuss the explicit forms of 
the evolution equations for the moments Mi , Mi , and S % 
and the wordlines z\ and z\, which depend only on these 
quantities and Q l K 

B. 1PN multipole and tidal moments 

In Sec. UH we defined the Newtonian mass multipole 
moments " M A (called simply M A there) as integrals 
over the body's mass distribution (|2.7[) . In so doing, 
we implicitly assumed that the Newtonian Poisson equa- 

2 Though the body-frame coordinates (sa_,v\) are only well- 
defined and physically meaningful in the buffer region B a, we 
will below make use of functions defined and evaluated at the 
spatial origin y A = 0. For this purpose, one can use any smooth 
extension of the coordinate system to y A = 0. While such ex- 



tion (I2.1|) was valid in all space, including the interior 
of the body. A similar approach can be taken at 1PN 
order, defining lPN-accurate multipole moments as inte- 
grals over a body's stress-energy distribution (as in (13. 12)) 
below), assuming that the 1PN field equations (|3.6I) are 
valid in all space. This was the approach taken in the 
original DSX formalism. 

As stressed by RF, one can also define a body's 1PN 
multipole moments without requiring the validity of the 
1PN field equations in the interior of the body. Instead, 
one need only impose the field equations in a vacuum 
buffer region Ba, a region of finite extent enclosed be- 
tween two coordinate spheres centered on the body; the 
moments can then be defined through the multipole ex- 
pansion of the 1PN metric in the region Ba- This allows 
one to consider objects with strong internal gravity, like 
neutron stars and black holes, as long as there exists a 
region Ba exterior to the object where gravity is weak 
enough for the 1PN field equations to be valid. 

Taking the latter approach, we assume the existence of 
a local coordinate system (s^,^), in the vicinity of the 
body A, having the following properties: (i) The range 
of the coordinates includes the product of the open ball 
\Ua\ < T%, for some finite radius ri, with an open inter- 
val of time (s\, s 2 A )- (ii) There exists a spatial region Wa 
(the worldtube) of the form ly^l < ri that contains all 
the body's stress-energy and/or regions of strong gravity, 
(iii) In the buffer region Ba {ri < \va\ < ri), 2 the coordi- 
nates (sa,Ha) are conformally Cartesian and harmonic, 
and the metric takes the 1PN form (|3.4[) . with potentials 
<S>a(sa, %)a) and Ca[ s a,Va) satisfying the 1PN vacuum 
field equations: 

V 2 $a = c~ 2 $ + 0{c- 4 ), 
V 2 ^ = 0(c- 2 ). (3.7) 

Under these assumptions, RF showed that the general 
solution for the potentials in Ba is of the form 

tended coordinates will have no direct physical significance, this 
is of no consequence to the validity of the formalism; see discus- 
sion in RF Sec. VC 



°° 1 f 1 

<S>a{sa,Va) = -Y, T A(-l) l Mk(sA)d L — + G L A ( S A)y 

£ Q ' k I S7-fi I 



(3.8a) 



^a{sa)Va + 



1 



;G%s A )y J i 



JjL 



Ca(. s a,Va) 



2(21 + 3) 

E T\ Z aM + YX l (s a ) yij + 0(c" 4 ), 



+ o( c - 4 ), 



(3.8b) 
(3.8c) 
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with 



j^Mf(s A ) - JL e *««S^>j {sA ) + |_l^<-^-i>(^) + {c-% 



l + l 



(sa) - -^fGr-'S^isA) + 0(c-*). 
I 



(3.9a) 
(3.9b) 



The potentials are parametrized by the following sets 
of (mulit-index) spatial tensors, which are STF on all 
their indices, and which are functions only of the time 
coordinate sa- First, M a (sa), with I > 0, are the 
body's mass multipole moments, which are defined with 
lPN-accuracy. Next are the current multipole moments, 
S a (sa), with I > 1, needed only to Newtonian accuracy. 
Together, the mass and current multipole moments con- 
tain all the information about the body's internal struc- 
ture that is encoded in the gravitational field it produces 
(at 1PN order). They are associated with the contribu- 
tions to the potentials that appear to diverge as \ua\ — >• 0, 
which can be referred to as the internal contributions. 
Also associated with such parts of the potentials are the 
internal gauge moments [i\ (/ > 0), so called because 
they contain no gauge-invariant information about the 
body. 

Associated with the parts of the the potentials that 
appear to diverge as \va\ — > oo (the external parts) are 
the tidal moments: the gravito-electric tidal moments, 
G\{s A ), with I > 0, are defined with 1PN accuracy (like 
M A ), and the gravito-magnetic tidal moments H A (s A ), 
with I > 1, arc defined with Newtonian accuracy (like 
S A ). The tidal moments contain information about grav- 
itational fields generated by external sources and about 
inertial effects associated with the motion of the local co- 
ordinate system. Finally, there are the tidal gauge mo- 
ments v\, defined for I > 1, 

The tensors Z a l (sa) and Y A (s A ) appearing in the 
gravito-magentic potential (|3.8cl) have been defined as 
useful shorthands for the expressions in (|3.9p . Unlike all 
the other moments just introduced, they are not STF on 
all their indices, but rather on all but the first. Eqs. (13. 9[) 
in fact represent their unique decompositions in terms of 
fully STF tensors; the 'inverse' relations are 



1 



-z 



jk<L- 



zf 



M 



iL 



l + l 



'<iL> 



and 



L _ yjk<L—l ai>jk 



m = Y- 



Gh 



= Y, 



<L> 



-Y 



30 L 



(3.10a) 
(3.10b) 
(3.10c) 

(3.11a) 
(3.11b) 

(3.11c) 



The relations (|3.10cj) and (|3. 1 lc|) arc implied by the har- 
monic gauge condition Q3.5p . 

In the case where the 1PN field equations (|3.6p are in 
fact valid in the interior of the body, the mass and current 
multipole moments can be defined by intergrals over the 
stress-e nerg y distribution in the volume of the body, as 
in DSX fill: 



Mi = 



j^yA^T" 



(3.12a) 



<L>rpjj 



2(21 + 3) 



(Z + 1)(2Z + 3) 



y n<L>fu 



|+0(c- 4 ) 



d 3 yA e jk<a ly L-l>j T tk + 0(c" 2 ). (3.12b) 



When considering a body with strong internal gravity, its 
interior cannot be modelled by a 1PN stress-energy distri- 
bution, and such integrals cannot be defined. Instead, we 
rely on the multipole exapansion of the potentials (|3.8[) 
in the buffer region B A to define the multipole moments 



Mi 



and S A , as well as the tidal moments G A and H A . 



Appendix E of RF [16| demonstrates the sufficiency of 
this method of definition by giving explicit formulae for 
the moments in terms of surface integrals of the poten- 
tials in Ba- 



Coordinate transformations and body-frame 
gauge conditions 



The 1PN metric Q3.4p harbors residual coordinate free- 
dom not fixed by the conformally Cartesian and harmonic 
gauge conditions. As a result, the mutlipole and tidal 
moments defined in the last section (not just the 'gauge 
moments' [i\ and is^, but rather all of the moments) are 
not unqiue and will vary with the choice of coordinates. 
To define a unique set of multipole moments for a given 
body, one must define a unique coordinate system that 
is specially adapted to the body. Thus we turn now to a 
discussion of 1PN coordinate transformations. 

In RF [r| , it was shown that the most general transfor- 
mation between two harmonic coordinate systems (s,y l ) 
and (t,x l ) in which the metric takes the 1PN form (|3.4p 
can be written as 
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x\s,y) = y i + z i (s) + 



^z kk (s)S' lj - a(s)S ij + e tjk R k (s) + ^z ij (s) 



-z i (s)S jk - z k (s)S ij 



t(s,y) = s + i [a(s) + i J (s)y J ] + i 
c c 



/?(s,y) + \a{s)y" + j^z j A (s)y jkk 



+ 0{c-% 



(3.13a) 
(3.13b) 



being parametrized by the following functions. The vec- 
tor z' l (s) provides a time-dependent translation between 
the spatial coordinates and is defined with 1PN accuracy. 
Each defined with Newtonian accuracy 3 are the rotation 
vector Ri(s), and the functions a(s) and j3{s,y) which 
transform the time coordinate. All of these may be arbi- 
trary functions of their arguments (within the bounds of 
their post-Newtonian scaling), except that f3(s,y) must 
be harmonic, V 2 /? = 0, in order to preserve the harmonic 
gauge condition. 

In the treatment of the iV-body problem, we will make 
use of one global coordinate system (t, x l ) and one body- 
adapted coordinate system (s A ,y A ) f° r eacn body A. 
The global coordinates (described further in the next sec- 
tion) are used to track the bulk motion of all the bodies, 
while the body-adapted coordinates are used in the local 
description of each body — in particular, to define their 
body-frame multipole and tidal moments. The transfor- 
mation between the {t,x l ) and (s Al y A ) coordinates will 
take the form (|3.13[) . with different 'worldline data' func- 
tions, 

V A = {z A (s A ),R A (s A ),a A (s A ),/3 A (s A , yA )} (3.14) 

for each body A. These functions may be viewed as con- 
figuration variables for the body-adapted frame, speci- 
fying its position, orientation, etc. relative to the global 
frame. 

In order to uniquely define the body-frame multipole 
and tidal moments, we must fix all gauge freedom in the 
body-adapted coordinates (s A ,y A ). It was shown in RF 
that this can be always be accomplished by imposing 
the following conditions, which define the body- adapted 
gauge: 

M A (s A ) = (3.15a) 
R A (s A ) = (3.15b) 
G A {s A ) = pl A {s a ) = (3.15c) 
ti(s A ) = v a (sa) = 0, I > 1 (3.15d) 

Eq. (|3.15a[) . setting the body- frame mass dipole M A to 
zero, fixes the body's center of mass-energy to the origin 



3 Though an 0(c~ 2 ) contribution to a(s) would contribute at 
0(e -4 ) in Eq. l|3.13bjl . this contribution can be absorbed into 
the function (3(s,y). 



of the spatial coordinates y A = 0. Setting the rotation 
vector R A to zero in Eq. (I3.15b[) fixes the orientation of 
the body-frame spatial axes to those of the global frame. 4 
If the extended body A were replaced by a freely falling 
observer at y\ = 0, Eq. (|3.15cl) would ensure that the 
time coordinate s A measures their proper time. Finally, 
the fact that all the gauge moments can always be set to 
zero by a coordinate transformation, as in (I3.15d[) . shows 
that they are pure gauge degrees of freedom. 

We can think of the body-adapted coordinates as defin- 
ing the body's local asymptotic rest frame [27], in which 
the effects of external gravitational fields and inertial 
effects have been removed as much as possible 4 . The 
body-frame moments — the mulitpole and tidal moments 
defined by (|3.8I) in the body-adapted coordinates — then 
take on the values that would be measured by a local 
comoving observer in the body frame. The body-frame 
multipole moments M A and S A are the quantities de- 
scribing the bodies' internal structure that will appear 
in the final form of the translational equation of motions 
for an iV-body system. 



D. The global frame 

To treat the orbital dynamics of a collection of sev- 
eral bodies A = 1 ... N, we consider N + 1 separate co- 
ordinate systems: one body-adapted coordinate system 
(s A , y\) for each body A, and one global coordinate sys- 
tem (t,x l ). We take these coordinate systems to have 
the following properties: (i) For each body A, the body- 
adapted coordinates (s A ,y A ) cover the worldtube W A 
and buffer region B A of the body and satisfy all the as- 
sumptions and gauge conditions outlined in Sees. 1111 Bl 
and MI CI (ii) The bodies' buffer regions B A are non- 
overlapping, (iii) The global coordinates (t,x l ) cover the 
buffer regions of all the bodies as well as the intervening 



4 In place of the condition d3.15bll . RF chose to set the gravito- 
magnetic dipole tidal moment H a (sa) to zero, which cancels 
leading-order Coriolis forces in the body-adapted frame and re- 
quires a non-zero value of the rotation vector R A (s A ). While 
this more completely effaces external gravitational and inertial 
effects in the body frame, Eq. (13 . 1 5bft leads to more simlifica- 
tions in calculations. The effects of these differing gauge choices 
cancel in all final results. 
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space; i.e. they cover the region B g — M\{J A W A , the the PN vacuum field equations f (13.7[) with A —> g). 
entire spacetime manifold M. except for the worldtubes. 

(iv) In the region B g , the coordinates are confor- The final assumption allows us to write down the fol- 

mally Cartesian and harmonic, and the metric takes the lowing multipole expansion of the global-frame potentials 

form (|3.4p . with potentials $ g (i, x) and Q(t, x) satisfying in B g : 



N 



(-1) 



= - EE 

A=l 1=0 

N oo , w 

= -EE 



1 



1 







A=l 1=0 



l\ g 



Z g L A (t)d L 



\x-z A (t)\ ' 2c 2 ^ 
1 



\x - z A (t)\ 



M^ A (t)d L \x - z A (t)\\ | + 0(c- 4 ), (3.16a) 
+ 0(c- 2 ), (z*£(f) = o), (3.16b) 



This expansion is analogous to that for the body-frame 
potentials (|3.8p but has several important differences. 
Firstly, the potentials are written as a sum of contribu- 
tions from each body A; this is justified by the linearity 
of the field equations (|3.6[) . Each such contribution is 
parametrized by the body's global-frame multipole mo- 
ments: the mass multipole moments M^ A (t) are fully 
STF and lPN-accurate, and the tensors Z l g L A {t) are STF 
on all but the first index and Newtonian-accurate. Both 
sets of tensors are defined for I > 0. As these global- 
frame moments will only appear in intermediate stages 
of our calculations, we will not bother decomposing the 
tensors Z^ L A in terms of fully STF current and gauge mo- 
ments as in the body-frame case (|3.9a|) . The condition 
Z gA = in (|3.16bl) is equivalent to setting the would-be 
global-frame gauge moments /Zg A to zero. 5 The global- 
frame moments M^ A (t) and Z^ L A (t) are distinct from 
(though related to) the corresponding body-frame mo- 
ments M*((s A ) and Z A L (s A ). 

A second important difference with the body frame 
case is that the multipole expansions appearing here are 

5 This condition, along with the fact that the global-frame poten- 
tials all vanish as \x\ — > oo, reduces the residual gauge freedom 
in the global-frame metric to the group of post-Galilean trans- 
formations (the post-Newtonian Poincare group) [T^l . which are 



centered not around the spatial origin x l = but around 
the worldlines x % = z A (t). One can check that the poten- 
tials as written here still satisfy the 1PN field equations 
for any choices of these worldlines. Below, we will iden- 
tify the z A with the bodies' center-of-mass worldlines, 
which appear as parameters in the transformations from 
body-adapted to global coordinates (|3.13p . 

Finally, one can note that we have included, in each 
body's contributions to the potentials, only internal 
pieces (which appear to diverge as \x — z A \ — > 0) and not 
tidal pieces (which would appear to diverge as \x — z A \ — > 
oo). This gauge choice makes the global- frame metric 
tend to the Minkowski metric as \x\ — > oo, thus elimi- 
nating any tidal or inertial forces on the A-body system 
as a whole. Each body will still experience local tidal 
fields, but they will arise from the contributions to the 
potentials generated by the other bodies. We can intro- 
duce a set of global-frame tidal moments for each body 
A by rewriting the global-frame potentials, in the body's 
buffer region B A , as 

the coordinate transformations given by A3.13I I with n z* = h % = 
Ri = p = and a = n z 2 /2. 



°° 1 f 1 



+ G L g At)[x-z A {t)] L 



(3.17a) 



5?* 



(-l) l M^ A (t)d L \x - z A (t)\ + -J—G^Mb - z A (t)]" L 



0{c 



oo 1 

En 



(-i)'2*Vt)0 



^ Ay "'" h \x-z A {t)\ E 



+ Y^ A [x-z A (t)Y 



+ 0(c-% 



(3.17b) 



Here, we have absorbed the contributions to the poten- tials from the other bodies B ^ A into tidal terms for 
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body A. This defines the global-frame tidal moments 
Gg A (f) and Y^ A {t). They can be expressed in terms of 
the global-frame multipole moments of the other bod- 
ies B 7^ A and the worldlines z\ of all the bodies A by 
equating the expressions for the potentials in (J37TTJ) with 
those in (|3.16p : these relations are given in [section]. 

Now, as the global coordnates x 1 - 1 — (t,x l ) and the 
body-adapted coordinates y A = (sa,u\) are related by 
the coordinate transformation (13.13[) . the metrics in the 
global and body frames must be related by the tensor 
transformation law: 



1 A 

7 \1V 



dx p dx 



(3.18) 



This requirement allows one to determine both the form 
of the coordinate transformation between the two coor- 
dinate systems (|3.13l) and the relationship between the 
global- and body-frame multipole and tidal moments. 
Making use of the form (I3.4[) for the metric in terms of the 
potentials (in both coordinate systems) and the expres- 
sions for the body-frame potentials (|3.8p and the global- 
frame potentials ([3~P7j) . as detailed in RF [lj| , Eq. (l3TTg)) 
yields expressions for the body-frame moments in terms 
of the global-frame moments and the worldline data (or 
the inverse relations): 



{M L A ,S L A ) 
(G A ,H A ) 



Pa. 



(AC 



L ryih \ 

g: A^g,A) 



(G 



g,A> z g,A 



These moment transformations are presented in full de- 
tail in [sees]. 

By combining the transformation formulae for the tidal 
moments with the body-frame gauge conditions (|3.15p . 
one can solve for and eliminate the worldline data func- 
tions c^sa) and /3a(sa, Da)- The only remaining piece 
of the worldline data T>a P-14[) is the translation vec- 
tor z a (sa)- Recall that the body-frame gauge condition 
M A = (|3.15aj) . setting the mass dipole to zero, fixes 



the body's center of mass-energy to the body-frame spa- 
tial origin y A = 0. Setting y A — in the coordinate 
transformation (|3.13[) and eliminating sa, we see that 
x l = z A (t) encodes the body's global- frame center-of- 
mass worldline, where 



40) = Za(sa) 



(3.19) 



is the quantity zf(sA) expressed as a function of t, with 
the function s A (t) found by setting y\ = in Eq. p,13bp 
(see Eq. (jBip and discussion thereabouts). The transla- 
tional equation of motion for a body A, discussed in the 
next subsection, can be written in the form of a second- 
order ODE for the global- frame CoM worldline z A (t). 



E. Single-body laws of motion and translational 
equations of motion 

The single-body laws of motion are constraints on the 
lowest-order multipole moments of any body which gov- 
ern the exchange of energy, momentum, and angular mo- 
mentum between the body and the gravitational field. 
The laws of motion at 1PN order were first found by 
DSX, who derived them by using covariant stress-energy 
conservation at 1PN order in the interior of the body. 
The same laws of motion were later rederived by RF 
by using the 2PN (next-to-next-to-leading order in 1/c 2 ) 
vacuum Einstein equation in a buffer region surrounding 
the body, thus extending their range of validity to include 
bodies with strong internal gravity. 

The laws of motion are written in terms of the body's 
multipole and tidal moments as defined by the expansion 
of the 1PN potentials Q3.8P and are valid in any coordi- 
nate system in which the spacetime metric takes the form 
given by (|3.4p and (13.81) — not just in body-adapted coor- 
dinates. The results are 



1 00 1 

^ = ""a E u [ {l + 1} " M * + * " G ^l + ° (c ~ 4) ' 
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(3.20a) 
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(Z + l)(2Z + 3) 
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(3.20b) 
(3.20c) 
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Eq. (|3.20a[) shows that the mass monopole Ma is con- 
served at Newtonian order (0(c°)), but not at 1PN or- 
der. As discussed further in Sec. IIIIFI below, Ma con- 
tains 0(c~ 2 ) contributions from the internal energy of 
the body, which can vary as tidal forces do work on the 
body. The law of motion (I3.20c|) for the the spin S A is 
the same Newtonian-order tidal torque formula found in 

Eq. mnp . 

The law of motion (|3.20bp for the mass dipole M A 
governs the evolution of the body's total linear momen- 
tum. The body's translational equation of motion can 
be derived by applying (|3.20b[) in the body frame, i.e. by 
applying it to the body-frame dipole moment, as follows. 

At Newtonian order, the 0(c°) part of M\ in (|3.20b|) 
gives the net force acting on the body, as M\ is the 
body's total momentum (c.f. (|3.12a[l ). Since the body- 
adapted coordinates are chosen to be mass-centered 
(M A = 0), this net force must vanish in the body frame. 
This apparent equilibrium in the body frame is achieved 
by the balancing of gravitational forces from the other 
bodies with inertial forces, which are due to the fact 
that the body frame is accelerating with respect to the 
(asymptotically) inertial global frame, along the world- 
line z A {t). Both of these effects are accounted for by the 
body-frame tidal moments G A ; from the O(c ) part of 
[ref], we have 

G l A — Gg,A ~ z\ + 0{c 2 ), 
G L A = Gl A + 0{c-% (l>2), 

B^Ak=0 1 1 

Using these relations in (I3.20b[) , the requirement of equi- 
librium in the body frame, M A = 0, determines the equa- 
tion of motion for the worldline z A (t): 

M ^t = E E ^ M tGf L A + 0(c" 2 ) (3.21) 

B^A 1=0 

00 

= E E ik M * M « d ^w^-\ +0{c ~ 2) - 

B=£Ak,l=0 ' ' 1 1 

This matches the Newtonian equation of motion found 
above in (12 . 16[) . 

At 1PN order, the procedure is essentially the same, 
but more involved. One begins by setting M A = in the 
body frame, with M A given by (|3.20bl) . To arrive at a 
suitable form for the final equation of motion, one must 
then rewrite the body-frame tidal moments G L A and H A 
of body A in terms of the body- frame multipole moments 
M B and S B of the other bodies B ^ A and the worldline 
data T>c for all the bodies C; the details of this procedure 
are presented in Appendix [B] 

In the end, one arrives at an expression for the accel- 
eration z\(t) of the lPN-accurate global- frame center- 
of-mass worldline z A (t), defined by (|3.19p . (As in the 



Newtonian case, the acceleration term, describing iner- 
tial forces in the body frame, emerges from the trans- 
formation laws for the body-frame tidal moments.) The 
expression depends only on the body-frame mass and cur- 
rent multipole moments M B (t) and S B (t), e the global- 
frame worldlines z l B (t), and the time derivatives of these 
quantities, for all bodies B: 

z A (t) = F A {z B ,z B , M&M&MZ, S%,S%). (3.22a) 

Similar (though simpler) manipulations applied to the 
laws of motion (|3.20ap and (|3.20c[) allow one to write 
equations of motion for the mass monopole M^ft) and 
spin S A (t) in terms of the same variables: 

M A (t) = F A {z B ,z B ,M%,M%), (3.22b) 
S A (t) = r A (z B ,M^). (3.22c) 

The explicit forms of the translational equations of mo- 
tion (|3.22ap and the mass and spin evolution equations 
(|3.22b|) and (|3.22c|) will be given for the M 1 -M 2 -S 2 -Q2 
case in Sec. IIII F[ and are given in the fully general case 
in RF (l|| as corrected by an upcoming erratum. 

To arrive at a closed set of evolution equations for 
the quanities z A (t), M A (t), and S A (t), for all bodies A, 
the equations of motion (|3.22[) must be supplimented by 
equations for the multipole moments M£(t) and S£(t) 
for I > 2. Finding such equations will require a model for 
the bodies' internal dynamics, which will be addressed in 
Sec. ED 



F. M1-M2-S2-Q2 truncation 

We have presented above the formalism for treating 
the 1PN dynamics of a collection of many bodies, each 
with arbitrarily high-order mutlipole moments. Here, we 
apply that formalism to the two-body system discussed 
in Sec. II Dl with a body 1 having only a mass monopole 
moment M\ , and a body 2 having a mass monopole M 2 , 
a current dipole, or spin, S\ = S 1 ', and a mass quadrupole 
M% = Q 2 = Q %J ■ More precisely, we truncate the in- 
ternal parts of body-frame multipole series (|3.8|) for each 
body according to 

*1M = -^ + 0(c- 4 ), 

m\ 



6 Here, and throughout, Mi(t) and S^(t) are the body- frame 
moments AI^(sa) and S a (sa) expressed as functions of t at 
y A = 0; c.f. Eq. (I B 1 b and surrounding discussion. 
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neglecting the moments My for I > 2 and Sy for I > 1, 
and 



1 

1? 



Q iJ '%|y 2 | + 0{c-% 



1 

I'd 



neglecting the moments M 2 for I > 3 and Sj" for I > 
2. (The external parts of these potentials will be just 
as in (|3.8p . with arbitrarily higher-order tidal moments.) 
These expressions for the body-frame potentials define 
the body-frame moments My(s\), M2(s2), Q l -'(s2), and 
S i (s2) as functions of the body- frame time coordinates s± 
and s 2 - These moments can be expressed as functions of 
the global time coordinate, written My(t), M^it), S l (t) 
and Q l 3(t), using the coordinate transformation (13. 13)) 
with y\ = (c.f. Eq. jBl)). For the bodies' global- 
frame CoM worldlines z\(t) and z 2 (t), we will use the 
definitions 



z l 2~z{, r = \z\, n l = z l /r, 



(3.23) 



as similar to (|2.17p . except that the worldlines are now 
defined with 1PN accuracy, and 

v{ = z{, w*=4 v^vi-vl (3.24) 

With these conventions in place, we can apply the 
laws of motion presented in Sec. IIII El to find the evolu- 
tion equations for the moments My(t), M 2 (i), and S l (t) 
and the global- frame center-of-mass worldlines z\ (t) and 
z\{€). The results involve only these quantities and the 
quadrupole Q % i(t). 

As body 1 has no higher-order multipole moments, the 
law of motion (I3.20aj) requires that its mass monopole My 
be constant in time: 

My = 0(c- 4 ). (3.25) 
The same law of motion applied to body 2 gives 

M 2 = -i (^Q ij G» - Q»G{\ + 0(c~% (3.26) 

where the body-frame gravito-electric tidal moment G\ 3 
is given as [rcf] by 



G^G^ 2 + 0(c- 2 ) 



Mr 



It is worth pausing here to compare this rate of change 
of the lPN-accurate mass monopole M2 with the rate 
of change Newtonian internal energy E 2 nt discussed in 
Sec. HTFl From (|2~33|t and ([3~2l)|) . we find that they are 
related by 



M 2 = or 2 (^ nt + 3t/ Q ) + 0(c" 4 ), (3.28) 
where Uq is the Newtonian gravitational potential energy 
associated with the quadrupole-tidal interaction given by 
(|2.28[) . The mass monopole M 2 thus contains contribu- 
tions not only from body 2's internal energy but also 
from the tidal part of its external gravitational potential 
energy. Though Newtonian internal energy is not a well- 
defined concept for strongly self-gravitating bodies, we 
will simply take the relation 

M 2 = "M 2 + c- 2 (E' 2 nt + 3U Q ) + 0(c- 4 ) (3.29) 

to define the quantity E 2 nt in the 1PN context, with "M2 
being the conserved Newtonian-order rest-mass contri- 
bution. This partitioning of M2, which is fully consis- 
tent with the equation of motion (|3 . 26[) , given (|2.34[) and 
(|2.28p . will be useful in our discussion of internal dynam- 
ics below. (A thorough discussion of the ambiguity in the 
total mass-energy of a body by an amount of the order 
of its tidal potential energy, and of tidal heating in GR, 
can be found in Ref. [IBJ].) 

The evolution equation for the spin (|3.20cD gives the 
following tidal torque on body 2: 



Jjk/ 



WG k 2 a + 0(c- 2 ), 



(3.30) 



with G| a being given by (|3.27|) . This coupling between 
the spin and the quadrupole, which is a purely Newtonian 
effect, is the essential reason we cannot (in general) ignore 
the spin-orbit coupling terms in the 1PN translational 
equations of motion for bodies with quadrupole moments. 

Finally, working from the laws of motion (|3.20b[) for 
the mass dipoles, we can apply the procedure outlined 
in Sec. ITnEl and Appendix [5] to the My-M 2 -S2-Q2 sys- 
tem to find the translational equations of motion for the 
worldlines z\ and z\. The results are 



M l z\{t)=Fl M + Fl s + Fi 



F, 



2,M 



■n <ij> +0(c- 2 ). (3.27) 



M 2 zm 

with the monopole contribtions 



F 2,S 



2,Qj 



(3.31a) 
(3.31b) 
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vi M 2 . 1 M 2 

F 1M = — n + — — 1 n 



c r 



2 ..2 3, a -, 2 5M X 4M 2 



2 V - v x - -(n w 2 ) 



Mi 



1 Mi 



2,M 



+ ^n a (4<-3« 2 Q )|+0(c- 4 ), 
+ u i n (4^-3<)} +0(c- 4 ), 



the spin contributions, 



F* tg = — -2 e abc S c S ai (4v b - 6n bd v d ) - 6n ai v b + 0(c~ 4 ), 



F. 



2.S 



1 M l abcgc 



38 ai (n M v d -u b )+6n a V + 0(c~ 4 ) 



(3.31c) 
(3.31d) 

(3.31e) 
(3.31f) 



and the quadrupole contributions, 



' c 2 I 2r 4 



-2n a c5 bi 



-Q°N 5n 



2 2 7, c c . 2 47Mx 24M 2 



2«<'-ttf--(n t# 



5r 5r 



1 2 V 27 2r r 



1 S 7lf r 
3M lAab( n abi + 2n *S U ) j + (9(c- 4 ), 



2.Q 



4r 2 

3Mi 
2r4 Q 

+2n a <5 bi 



(3.31g) 



1M 


- 5n abl 







o 2 2 7, c c , 2 8Mi 6M 2 



r 2r 



r r 

i ah i r~,<3C2/r> b c bc\ 

+ n v + on (Zv v 1 — v 2 ) 



+v l (5n abc - 2n a 5 bc ){^v c 2 - 3t$ + n% 2 (w 2 - 2uJ) + 5 bl n c [(5t£ - 4t>?)«§ - 6u a <] 



3M! - at 



v"(2n M - S al ) + 8 m n bc v c - 2n a V + 0(c" 4 ) 



(3.31h) 



(It should be noted that occurences of 5* in the equa- 
tions of motion have been replaced by (|3.30j) and in- 
cluded in the quadrupole contributions.) The monopole 
contributions (|3.31cl3.31d[) give the well-known Lorentz- 
Droste-Einstein-Infeld-Hoffmann accelerations, and the 
spin contributions (I3.31cl3.311]) give the well-known 1PN 
spin-orbit terms [151 ] . The quadrupole contributions 
p.31g|3.31h[) have been derived previously by Xu, Wu, 
and Schafer [17j . though our results disagree with theirs 
in several terms; we have not been able to pin down the 
source of the disagreement, but note that their monopole 
contributions also disagree with the standard EIH result 
in a way that cannot be an artifact of differing gauge 
choices. Our results also disagree with the final results 
of RF [lH, but agree with their corrected results given 
in an upcoming erratum. The strongest indication of the 
correctness of our expressions for the EoMs is the fact 
that, unlike the results in [l6|, [I?]], they are consistent 
with the conservation of the binary system's total linear 



momentum, as discussed in Sec. IIV Cl below. 

IV. SYSTEM MULTIPOLE MOMENTS AND 
CONSERVATION LAWS 

In Sec. IIIIB1 we defined the multipole moments of a 
single body through the multipole expansion of the met- 
ric in a vacuum buffer region surrounding the body. The 
same procedure can be applied to a collection of several 
bodies to define multipole moments for the entire sys- 
tem. Applying the general laws of motion discussed in 
Sec. MI El to these system multipole moments will allow 
us to formulate conservation laws for the energy, momen- 
tum, and angular momentum of an isolated TV-body sys- 
tem, expressed as constraints on the worldlines and mul- 
tipole moments of the constituent bodies. These conser- 
vation laws can serve both as a consistency check for the 
equations of motion given in Sec. lIIIFl and as a means to 
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specialize the equations of motion to the system's center- 
of-mass frame. 

A. General formulae 

We have already discussed, in Sec. IIII D[ a form for 
the metric generated by a system of N bodies. Using 



the global coordinate system (t, x l ), we expressed the po- 
tentials parametrizing the metric as a sum of multipole 
expansions for each body A, written in terms of the bod- 
ies' global-frame multipole moments M^ A and Z^ L A and 
Newtonian-order worldlincs z\: 



(-1)' 



EE „ 

A 1=0 
oo ( 

EE 



A i=o 



11 



Z l L A d L- 



1 1 

x — za | 2c 2 

1 



dl (M^ A d L \x - z A \) 



+ 0(c- 4 ), 



x - Z A 



+ 0(c-% 



(4.1a) 
(4.1b) 



with Z^ A = being enforced as a gauge condition 
(c.f. (|3.16p ). This solution for the metric was constructed 
to be valid in the region B g , which extends out to spatial 
infinity. 

In a region far outside the system, we can rewrite these 
expressions for the global-frame potentials to mirror the 



forms p.8p used to define the mulitpole moments of a 
single body, with multipole expansions about the global- 
frame origin x l = 0: 



E 

1=0 



(-1) 



H 1 Ivl sys u L |^| T ^ 2 



1 1 

\x\ 



(Z + 1)(2Z + 3)^ sys u \x\ 2 



+ 0{c-% 



00 r_iv i 



1=0 



11 sys "\x\ 



(4.2a) 
(4.2b) 



with 



yiL — j .lift _ ]i<a, qL~l>j , ^ \ X i .,, ,,L i 

sys l + 1 sys l ^ ^sys ^ 2 l + 1 



F«»^-i>+0(c-*). 



(4.3) 



These expansions define the system multipole moments 
M S y S and S^ ys and the gauge moments /x^ . The tidal 
terms present in (|3.8p are absent here, as the global-frame 
potentials vanish as \x\ —> oo (c.f. Eq. (|4.1jl ). To compare 
the metric (|4.2p here to the metric (|4.ip above, we must 
express them in the same gauge. We have chosen the 
gauge that enforces Z^ A = in (|4.ip . which will result 
in nonzero values for the system gauge moments p^ ys in 
(EI). 



Since the potentials given by (|4.ip and by (|4.2p rep- 
resent the same metric in the same gauge, they should 
be explicitly equal. This condition will allow us to solve 
for the system multipole moments appearing in (|4.2p in 
terms of the individual bodies' global-frame multipole 
moments and worldlines appearing in (|4.ip . 



Considering first the vector potential Q, we can use 
the Taylor series 



k=0 



to rewrite (|4.1bp in the form 

°° (-l) l+k 1 

Q = -E E —^^IdLKT^ 
A Lk=0 1 1 



EEE 

A p=0 k=0 



z 



<P-K K> 



(4.4) 



p\ k\(p-k)\ s 



]A z A a P 



In the second line, we have relabeled the multi-indices 
according to LK — > P, adjusted the summations accord- 
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ingly, and used the fact that dp\x\ 1 is STF. Renaming 
P — > L and comparing this with (I4.2b[) gives an expres- 



sion for the tensors Z. 



iL . 



71 L 



EE 

A k=0 



V. 



{l-k)\k\ 



7 KL-K nK> 
^ g ,A Z A ■ 



(4.5a) 



The current moments and gauge moments can then be 
found from formulae analogous to (|3.10al) and (I3.10b|) : 



9^ = 

sys 
Msys = 



_yjk<L—\ ai>kj 
A sys 

Z jjL 
sys 



(4.5b) 
(4.5c) 



The scalar potential $ g can be manipulated in a similar 
manner, using the Taylor series (|4.4I) . to find formulae for 
the system mass multipole moments M^ ys . The details of 
this more involved procedure are given in Appendix [C] 
The result is 



M L 

sys 



EE 

A k=0 



kl(l - k)\ 



M <L-K K> 
1VI g,A A 



1 1 «2 ( 9M ]<L~K K> 3 M <L-K K>jj 



c 2 2(2/ + 3) 
1 21 + 1 
(Z + l)(2Z + 3) 



/ i s i vs + 0(c- 4 ). (4.6) 



The gauge moments fi^ ys appearing here can be found 



from (|4~5|) . 

In summary Eqs. (|4.5p and (|4.6[) give the total multi- 
pole moments Af s ^ s and S^ ys of an A^-body system, de- 
fined in the global frame (i, x l ), in terms of the individual 
bodies' global-frame multipole moments M£ A and Zl L A 
and worldlines z l A . In the following subsections, consid- 
ering the two-body M1-M2-S2-Q2 case, we will use these 
results along with the moment transformation formulae 
from Appendix [B] to write the system's mass monopole 
M sys and mass dipole M s l ys in terms of the body-frame 
moments (Mi, M 2 , S l , Q y ) and the worldlines z\ and z l 2 . 
We note that a similar procedure can be applied to find 
S l sys = ^ k {M x z{v\ + M 2 z J 2 v$) + Si for the system's to- 
tal (Newtonian) angular momentum. The system's 1PN 
accurate mass quadrupole, which will be needed for the 
calculation of the gravitational wave signal from the bi- 
nary system, can also be calculated from Eq. 



B. System mass monopole 

Specializing the general formula (|4.6p for the system 
mass multipoles to the monopole (I = 0) case, and using 



the M1-M2-S2-Q2 truncation, we find the binary system's 
total lPN-accurate mass monopole to be 



M g ,i + 
1 d 2 



M„ 



1 



1 



;A*sys 



0(c~ 4 ). 



, 6dt 2 (M gil z 2 + M g , 2 z 2 

Using (|4.5[) for fi sys and the formulae in Appendix IB II 
relating the global- and body-frame multipole moments, 
we can rewrite this expression in terms of the body-frame 
multipole moments and the CoM worldlines: 



Ms- 



M 1 + M 2 + — 



1 (M x v\ 



>-Z u 2 



M X M 2 



2 2 
2U Q )+0(c- i ), 



where the tidal potential energy Uq, as in (|2.28[) . is 



3Ml 



If we rewrite the mass monopole of body 2 as M2 = 
"M 2 + c- 2 (E 2 nt + 3U Q ), as in (j3~2^1) . wc^find that the 
1PN contribution to the system mass monopole is exactly 
the system's total Newtonian energy E given by (|2.32l) : 



M sys = Ah 



"Af 2 + c- 2 E + 0(c~ 4 ) 
2 



(4.7) 



This is a further validation of the decomposition of the 
total mass monopole M 2 in Eq. (|3.29p . The constancy of 
M sys , required by the law of motion (|3.20aj) as applied 
to the entire system in the global frame (for which there 
are no tidal moments), then follows from the constancy 
of E. 



C. System mass dipole 



Taking the I = 1 case in the general formula 
the system's lPN-accurate mass dipole: 



gives 



3 

lo' 



10 



df l2M»j4 + M gtl z?>+M s , 2 z? 



+ 0(c~ 4 ) 



Using (|4.5I) for fi z sys , (|3.30l) to replace an occurence of 

S l , (|3T2T))) to replace M 2 , (ET2"g)l for Uq, and the moment 
transformations from Appendix IB II this becomes 
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JrJ -sys 



Mi z 



M 2 z\ + \ 



M x vl M X M 2 Uq 
~~ 2 2~T~ + ~Y 



+4 



M 1 M 2 
2r 



Uq 
2 



E. 



int. 



3Mi 

~2~^ 



Q 



+ 0(c- 4 ) 



(4.8) 



From the law of motion (I3.20b[) as applied to the en- 
tire system in the global frame, for which there are no 
tidal moments, we see that M s l ys should vanish; this is 
a statement of total momentum conservation. By differ- 
entiating (I4.8[) . order reducing as appropriate, using the 
full 1PN translational equations of motion p.31[) in the 
Newtonian terms and their Newtonian parts in the 1PN 
terms, and also using (|2~34|) and ([3T30]) for E 2 ut and Si, 
we find that indeed M s ' ys = 0. This is an important check 
of the correctness of our expressions for the equations of 
motion and of the consistency of the formalism (which is 
not satisfied by the EoMs given in Ref. [H, Il7j). 



V. ORBITAL DYANMICS IN THE SYSTEM'S 
CENTER-OF-MASS FRAME 

A. Equation of motion of the relative position 



The conservation of momentum allows us to reduce the 
problem of solving for the two worldlines z\ (t) and z\ (t) 
to solving for just their separation z l (t) = z\(t) — z\(t) 
in the binary system's center-of-mass (CoM) frame. We 
can define the CoM frame to be that in which the 1PN- 
accurate mass dipole vanishes, 



with 

V = 
V, = 



, u 1 M 3 
V(X2 -Xi) "T ~ o ~a 3 

2r 2 ^ '" J ' M 



Q v nj + ^e ijk v j S k , 



Xl Tjiint 



(5.4) 



and with the new notation 

M = Mi + "M 2 , xi = Mi/M, X 2 = n M 2 /M, 

fi = Mi °M 2 /M, n = X1X2 = n/M. (5.5) 
To find the acceleration of the relative position in the 
CoM frame, we can simply subtract our above results 
p.3ip for the individual accelerations: 

a 1 = z l =z\ — z\. (5.6) 

The resulting expression depends only on z 1 , v\, v\, M\, 
M 2 , Q ij , S\ and Ef\ As v\ and v\ appear only in 1PN 
terms, we can replace them with their Newtonian values 
in the CoM frame, 

v[ = - X 2v l + 0(c- 2 ), vi = X iV i + 0{c- 2 ), (5.7) 



Ml ys {t) = 0, 



(5.1) 



so that the system's center-of-mass(-energy) is at rest at 
the global-frame spatial origin. This fixes all remaining 
(post-Galilean) coordinate freedom in the global-frame 
metric. Using Eq. (14.8[) . this condition can be used to 
solve for the worldlines z\ and z\ in the global CoM frame 
in terms of the relative position z l (working perturba- 
tively in c~ 2 ); one finds 



-X2Z l + c- 2 (Vz l ~V l ) +0(c" 4 ), 
0(c- 4 ), 



z\ = xiz l + c- 2 (Vi 



V 1 



(5.2) 
(5.3) 



from differentiating the 0(c ) parts of ()5.2|5.3j) . After 
defining one last shorthand, 

r = n a v a , (5.8) 

we can write our result for the lPN-accurate CoM-frame 
relative acceleration as follows: 



a 1 = a l M + a l s + a' Ql 



with the monopole contribution, 



(5.9a) 



M 



1 M 



TI 71 ~ ~^2 ~2 ^ " 



2 3/7 . 2 K M 



(l + 3r?K - -^-r 2 -2(2 + 77) 

2 r 



2(2 - 77)7V \ + 0(c~ 4 ), 



(5.9b) 



the spin contribution, 



e abc S c 

C 2 X2T 3 



[(3 + X 2)v a 8 u - 3(1 + X2)rn a 5 bl + 2n al v b ] + 0( C - 4 ), 



(5.9c) 
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and the quadrupole contribution, 
3Q ab 



[5n abi - 2n a S bt ] + 



i Q 



ab 



n abi ( B lV z + B 2 r z + B 3 



M 



M 



-n a S bl [B^v 2 + B 5 r 2 +B e —j+ B-rn ab v l + B 8 n a v bl + B g rn al v b + B w v ab n l + B n rv a S bl 
[Bi 2 n ab v i + B 13 rn aM + B lin m v b + B 15 v a S bl + B 16 rn a S bl ] + %- [B l7 n aU + B ls n a S M ] 



rpint 

R 2 
-£>- 



>19" 



n' +0(0: 



(5.9d) 



with coefficients 



B, = ~(1 + 3 V ), B 2 = B 3 = H( 5 - 2*2), s 4 = A( 2 + 2%2 - 3*1), 

2X2 4 X2 X2 



15 
'2x~ 2 



^ : ^(2 -X2- xll B 6 = - A ( 8 - X2 - 3x1), S 7 = — (2 - 77), B, 

X2 



15 
— ( 

X2 



3 
2^ 



s -7-(7-2 X2 +3x^), 



= — 1 + X2 ), Bio = 7; — , Bn = - — 5 - 4x2 - X 2 )> B 12 = -- — (4 - X2), B u = — , 

2X2 2x2 2x2 ~ 2x2 2 



B14, = — , Bi 5 — — , B 16 — — (1 - 2x2 - xl)) Bn = -, Bi S — ' Big — 1. 

X2 X2 X2 4 2 



(5.9e) 



In this form for the CoM-frame orbital EoM, we 
have used (I3.29P to write the total lPN-accurate mass 
monopole M 2 in terms of the (constant) Newtonian mass 
" M 2 , the internal energy E 2 nt , and the tidal potential 
energy Uq (giving a contribution to B 3 ). This decompo- 
sition is useful in formulating an action principle for the 
orbital dynamics (as in the next subsection), as E 2 nt is 
independent of the orbital degrees of freedom, while M 2 
is not. 



B. Generalized Lagrangian for the orbital dynamics 



The monopole contributions (|5.9bp to the 1PN CoM- 
frame orbital EoMs are known to be derivable from the 
Lagrangian 



-M 



fJ,V 



r 



1 - 377 4 
— - — v 



M . 

-_((3 + „),r 



• 2 M 

rjr 

r 



(5.10a) 
+ 0(c- 4 ), 



(see e.g. [28[). The spin contributions (|5.9cp can also be 
derived from an action principle, but with a generalized 
Lagrangian (one depending not only on the relative posi- 
tion z % and velocity v % = z\ but also on the acceleration 



a 1 = z 1 ) given by adding 



JLs — O V 



2M 



Xi 



+ 0(c- 4 ) (5.10b) 



to (I5.10a|) (see e.g. [29|). Applying the generalized Euler- 
Lagrangian equation, 



d d d d 2 d 
dz l dt dv l dt 2 da 1 



C = 0, 



(5.10c) 



to C = Cm + Cs, and using Mi = M 2 = S l = (which 
replaces (|3. 2513.2613. 30p in the case with no quadrupole), 



one recovers the EoM a 1 = a\ 



M 



from (|5.9bl5.9cl) 



We have found that the quadrupole contributions to 
the orbital EoM (|5.9dP can also be encoded in a general- 
ized Lagrangian. To determine the necessary additions to 
the Lagrangian, one can proceed by guesswork, using the 
known Newtonian Lagrangian (|2.27p . and writing down 
all possible lPN-order scalars that can be formed from 
the relative position z % and velocity v l , the total (Newto- 
nian) mass M, and linear factors of the quadrupole Q % \ 
its time derivative Q 1 ^ , and the internal energy E 2 nt ; in- 
cluding dimensionless coefficients Ai~Ag for each such 
term, we have 
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r - 3Ml o^<> + 1 J M n« b 



n Q6 ( An; 2 + A 2 t 2 + A 3 — I + A 4 v a& + A 5 fn a v 



/*\Q>b [a „a b 1 a '„a61 , rni 



tint 



A 8 v 2 + Ag- 



AI 



+ 0(c" 4 ). (5.10d) 



Since the spin-orbit terms (|5.10bl) require the accelera- 
tion a 1 , one might expect that terms with factors of a' 
and also Q lJ should be included here; we find, however, 
that such terms are not necessary. The only further term 
allowed by general considerations but not included here 
is E 2 nt f 2 , as recovering the EoM (|5.9|) requires its coeffi- 
cient to be zero. 

By applying the Euler-Lagrange equation (|5.f Ocj) to 
the generalized Lagrangian C = Cm + Cs + Cq , using the 
evolution equations (13.30)) and (|2.34[) for time derivatives 
of S l and E^, one finds an EoM of the same form as (|5.9[) 
but with coefficients Bi~Big in (|5.9d[) given as functions 
of the Lagrangian coefficients A\-Ag and the mass ratios 
Xi and X2- Setting these coefficients equal to the values 
for Bi-Big given (|5.9el) gives a system of 19 equations for 
the 9 unknowns Ai~Ag, which has the unique solution 



A 1 
A 3 
A 5 
A 7 



3xi 



:i+3xi) 



-^(3 + X 2 ) 
3r/ 

'T' 



At = 5*1, 
2 ' 

b 2 



As 



X? 



Xi- 



(5.10e) 



Thus, the action principle (I5.10|) reproduces the 1PN 
CoM-frame equation of motion (|5.9I) for the relative po- 
sition — if we also make use of the evolution equations 
(1X30)) and (|234| for the spin S i and internal energy E' mt 
of body 2. In the next section, we discuss an action prin- 
ciple that leads to a closed set of evolution equations for 
the binary system in the adiabatic approximation. 



VI. INTERNAL DYNAMICS IN THE 
ADIABATIC APPROXIMATION 

A. Euler-Lagrange equation for Q ab 

We have just seen that the CoM-frame orbital EoM 
(|5.9[) , the evolution equation for the binary's 1PN- 
accurate relative position z l (t) — z\{t) — z\(t), can be 
derived from the action principle (I5.10[) . We now seek to 
extend this action principle to encorporate the internal 
dynamics of the deformable body 2 in the case where the 
quadrupole moment is adiabatically induced by the tidal 
field. 



In Sec. Ill G[ we saw how the adiabatic evolution of the 
quadrupole can encoded in a Newtonian action principle; 
varying the action (12.36)) with respect to the quadrupole 
Q l i gives 



0(c 



\^± n <v> + (c 2 ), 



(6.1) 



for the Newtonian-order quadrupole (c.f. (|2.35)) V We 
also saw that, with the quadrupole given by (|6.1|) . 
the spin evolution equation becomes S l = 0(c~ 2 ) 
(c.f. Eq. I2.46[) . so that body 2 experiences no tidal 
torques. For this reason, in the adiabatic case (unlike 
in the general case), we can specialize our analysis to 
the case of zero spin without generating inconsistencies, 
which we will do for the remainder of this section. 

We have found that the simple Newtonian Lagrangian 
(|2.36j) can be extended to govern the lPN-accurate adia- 
batic evolution of Q lJ in a relatively straightforward man- 
ner. We consider the following Lagrangian: 



C 

^orb 



= £ orb + 4 nt + o( c - 4 ), 



C 



C 



int. 



M 



4A 



U a0 Q 



abr^ab 



V 



ab/\ab 



WE" 



Q Q 



(6.2a) 
(6.2b) 



by (j5TUa|) and with U ab (z,v), 
being the coefficients appearing 



with Cm{z,v) given 
V ab (z, v), and W(z,v 
in Cq (|5. 10d[) . We have postulated (motivated by sym- 
metry considerations) that the internal Lagrangian C 2 nt 
can still be taken as a simple quadratic in Q ab (|6.2bj) . gen- 
eralizing the Newtonian Lagrangian (|2.36[) only by using 
the lPN-accurate value of Q ab in place of its Newtonian 
value and by using the lPN-accurate value for the tidal 
deformability A. The internal energy E mt appearing in 
C OI h (|6.2ap . which is needed only to Newtonian order, 
will still be given by 



E int = -±-Q ab Q ab + 0(c- 2 ), 
4A 



(6.3) 



up to a constant, as in (I2.40[) . To avoid explicitly intro- 
ducing an additional constant contribution to the inter- 
nal energy, we can absorb any such contribution into the 
constant 'Newtonian' mass monopole n M 2 (c.f. Eq. l3.29[) . 

Treating the relative position z l {t) and the quadrupole 
Q ab (t) as the independent dynamical variables in the La- 
grangian (|6.2j) . we still recover the orbital EoM (|5.9p from 
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the Euler-Lagrangian equation for z % (with S a = 0), and 
that for Q ab , 



8 



d d 



dQ ab dt dQ al 



C = 0, 



gives 



\ab 



2A(1 + W){-U ab + V ab ) + 0(c~ 4 ) 



(6.4) 



■5) 



Using the coefficients U ab , V ab , and W from the La- 
grangian (I6.2al5. 10|) . we find that 



Q ah = XGf + 0(c" 4 ) 



3.6a) 



where is the body-frame gravito-electric tidal mo- 
ment defined in Sec. IIIIBI and calculated in Appendix 

m 



G ab = ^XlM_ n< ab> 



1 3 X iM 



2 5x| .2 5 + X i M\ 

2~ r 2 V ' " '"' 



+V<ab> ~ (3 - X2)™<a' u b> 



■0(c" 



(6.6b) 



The relation (|6.6a|) between the body-frame tidal mo- 
ment and the adibatically induced quadrupole Q , 
derived here from the Lagrangian (|6.2j) . is the same rela- 
tion typically used to define the adiabtic approximation 
(e.g. in [2l|). The explicit expression for the tidal mo- 
ment in (|6.6b[) matches that given in Ref. [3(| when the 
latter is specialized to the CoM frame via (I5.7[) . 



and with the dimensionless coefficients 

O = (l-3r?)/8, 

0! = (3 + 7?)/2, 

02 = fi/2, 

03 = -1/2, 

6 = ( X i/2)(5 + X2 ), 

6 = -3(l-6x 2 +xD, 

6 = -7 + 5x2- 



(6.9) 



While this form for the Lagrangian has been derived in 
harmonic gauge, we note that (some) other gauge choices 
lead to a Lagrangian with the same terms as in (16.71) but 
with different values of the and £ coefficients. In par- 
ticular, the Lagrangian obtained (via a Legendre trans- 
formation) from the EOB Hamiltonian including 1PN 
tidal effects proposed by Damour and Nagar [2l| has this 
form, except that their work has not provided a value for 
the coefficient 6- In Appendix [S] we derive the canoni- 
cal transformation relating the EOB Hamiltonian to the 
harmonic-gauge Hamiltonian, which fixes the value of 6 
in the EOB Hamiltonian. 

The orbital EoM resulting from the Lagranigan (|6.7[) 
is given by 
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(6.10) 



+ fV 
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with coefficients 



B. Reduced Lagrangian, equations of motion, and 
conserved energy 



By substituing the solution (|6.6p for the quadrupole 
into the Lagrangian (|6 . 2[) . we find a reduced Lagrangian 
for the orbital dynamics involving only the CoM-frame 
orbital separation z z (t): 
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1 + ^5 
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(6.7) 
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C3 


= 120!+ 


120 2 + 26 + 26 - 76 
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= 12(40o 


+ 6), 



(6.11) 



for general values of the Lagrangian coefficients, and with 

01 = -1 - 3ry, 

02 - 377/2, 

03 = 2(2 + 7?), 

04 = 2(2-7?), 

Ci - -3(2-x 2 )(l + 6x 2 ), 

6 = 24(1- 6x2 + X 2 ), 

Ca = 66 + 9x2-19x2, 

C4 = 6(2-x 2 )(3-2x 2 ), (6.12) 



2G 



in harmonic gauge. 

Finally, from the Lagrangian 
the conserved energy, 



, we can construct 



E = v l dC/dv l -C 
pv 2 \iM 
~2 V 
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1 + T 



fj, 



39 v 4 
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V 2 I (7, 
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A 



£3 
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,(6.13) 



which is a constant of motion of the orbital EoM (|6.10l) . 

As an application of these results, we can compute the 
gauge- invariant energy-frequency relationship for circular 
orbits. Using the relations r = 0, v 2 = r 2 w 2 , and a 1 = 
—ruj 2 n l for a circular orbit, the orbital EoM (|6.10l) can 
be solved perturbatively, working to linear order both 
in the post- Newtonian parameter 1/c 2 and in the tidal 
deformability parameter A (|6.8p . to find the radius r as 
a function of the orbital frequency uj. Combining this 
result with a similar treatment of the energy (I6.13[) . we 
can eliminate r to find E(u): 



E(u) = p(Muj) 2/3 



1 

~ 2 

+ /. 



M- 



Jm 



3Aw^3 
f M 5 / 3 
( Mm) 2 / 3 

9 + rj 



Iq 



Alj 4 
M~c 2 



~(e + e 1 + e 3 ) = . >4 . 

y(80 O + 201 + &) 

^(3 + 2 X2 + 3 X ^). 
6 



3.14) 



While the 9 and £ coefficients may take different values 
in different gauges, their combinations appearing here 
must be gauge- invariant. As E(ui) is independent of £2, 
we can check this result against those obtained from the 
EOB Hamiltonian of Damour and Nagar (see Appendix 
IX)) . and we find that they agree. 



VII. CONCLUSION 

We have derived the post-l-Newtonian orbital equa- 
tions of motion for binary systems of bodies with spins 
and mass quadrupole moments, at linear order in the 
spin and quadrupole, and shown that they conserve the 
total linear momentum of the binary. After specializing 
these results to the binary's center-of-mass-energy frame, 
we have found an action principle from which the orbital 
equations of motion can be derived. Finally, we consid- 
ered the case in which the quadrupole moment is adia- 
batically induced by the tidal field, giving a simplied La- 
grangian and equation of motion for this case, as well as 
the conserved energy function and the energy-frequency 



relationship for circular orbits. These results will be use- 
ful in calculating tidal effects in the GW signal from in- 
spiralling neutron star binaries, which will be addressed 
in future work. 
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Appendix A: Hamiltonian for the adiabatic orbital 
dynamics, canonical transformations, and 
comparison with EOB results 



As in Sec. lVI Bl we consider a Lagrangian for the CoM- 
frame orbital separation z l = z\ — z\ of the form 

„2 1 / A - 



£ 



V" 1 , , 

— + - I 1 + 

r 



6 



A 



0(c- 4 ), 



— 02 



(Al) 



in units that set M = 1, and with the Lagrangian having 
been rescaled by the symmetric mass ratio, C = C/rj. We 
have shown that the harmonic-gauge values of the 6 and 
£ coefficients are given by Eq. (16. 9|) . 

The (rescaled) momentum canonalically conjugate to 
z l is 



P 



dC 

dv l 



1 



2v i ( A 
40 O «V + — U +£1 — 
r \ r° 



2rv l 



6 



A 



0(c- 4 ) 



±5p l (z,v) + 0(c- 4 ). 
& 



(A2) 



From a Legrendre transformation of the Lagranigan 
C(z, v), we can construct the Hamiltonian: 



H(z,p) — p J v 3 



C(z,v) 



A*,p)+^|§^+o( c - 4 ) 



£(z,p) + 0(c- 4 ), 



(A3) 



having used 5p l (z,v) = Sp t (z,p) + 0(c 4 ). This gives 



v 1 ( A 
H = — -- 1 + ^r 



0o^ 4 + — 
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+ 0(c- 4 ) 
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(A4) 
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In Ref. [21|, Damour and Nagar consider an EOB 
Hamiltonian of the form 



1 

Heob = — 
V 



1 + 2ri(H eB - 1) 



with 



H, 



off 



1/2 



1 1/2 



(A5) 



(A6) 



Here, p r and p^ are the momenta conjugate to the po- 
lar coordinates (r, 4>) in the plane of motion, related to 
the Cartesian momenta used above by p r — n ■ p and 
p^/r 2 = p 2 —p 2 . The functions A(r) and B{r) are coeffi- 
cients in the EOB effective metric [3l|; A(r) completely 
encodes the energetics of circular orbits, while B(r) has 
effects only when p r ^ 0. Damour and Nagar proposed 
to incorporate Newtonian and 1PN tidal effects into this 
EOB Hamiltonian by adding tidal terms to the radial 
potential A: 

2A 



Mr) = l- ^--4^( 1 + ^ L )+ °( c ~ 6 )' ( A7 ) 
c r c r \ c r' 



for the quadrupole I = 2 case. They have computed the 
1PN tidal coefficient to be 



Oil = 2*2- 



(A8) 



While they did not propose to modify the potential B(r) 



from its point-particle value of B 
we find that such a modification, 



1 + 2/c 2 r + 0(c 4 



B(r) 
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1 



a4 



0{c 



(A9) 



for some coefficient (3q, is necessary to match our results. 
Expanding the EOB Hamiltonian with these values for 
the potentials, we find a Hamiltonian of the form (|A4[) 
with coefficients 



9i 

h 



(l + »7)/8, 
1, 

/3b, 

l + 'q + ax, 



(A10) 



instead of the harmonic-gauge coefficients in Eq. (I6.9[) . 

Without tidal effects, the EOB Hamiltonian and the 
harmonic-gauge Hamiltonian (which coincides with the 
ADM Hamiltonian at 1PN order) are known to be related 
by a canonical transformation |3l| . Considering a 1PN- 
order canonical transformation with generating function 
G, 



1 d_ 

dp' 



G(z,p) 



c 2 dz' 



G(z,p), 
(All) 



the Hamiltonian changes by 



1 



H^H+-{H,G} + 0(cr 4 ), 



(A12) 



with only the Newtonian part of the Hamiltonian con- 
tributing in the Poisson bracket. We find that the most 
general generating function G that preserves the form of 
the Hamiltonian including tidal effects (|A4I) , changing its 
coefficients but adding no new terms, is of the form 



G = (z-p)[ 7ip- 



1 



A 



72- +73^ 



(A13) 



with arbitrary constant 7 coefficients. The changes in the 
Hamiltonian coefficients induced by the canonical trans- 
formation are 



A6 
A6» 2 

A?3 



-7i) 
7i - 72, 
27i + 72, 

72, 

671 - 73, 

1271 + 673, 

672 + 73, 



(A14) 



If we set the EOB Hamiltonian coefficients (lAlpp equal 
to the harmonic Hamiltonian coefficients (I6.9[) plus the 
transformation parameters (|A14[) , we find that this (re- 
dundant) system of equations has a unqiue solution. The 
coefficients in the canonical transformation (|A13|) must 
be 



71 = — *?/2, 

72 = (2 + t?)/2, 

73 = 2-9 X2 /2 + 2 X ^ 



(A15) 



with the values for 71 and 72 matching those computed 
in Ref. [3l|, and the parameters in the EOB potentials 
must be 



"i = 5x2/2, 
A, = 3(3-5?/). 



(A16) 



The value for ot\ matches that given by Damour and Na- 
gar, and the value for /3q can be used to extend the range 
of validity of their EOB Hamiltonian to non-circular or- 
bits. 



Appendix B: Moment transformations and 
translational equations of motion 

We present here the formulae that relate the body- 
frame multipolc and tidal moments and the global-frame 
multipole and tidal moments. They are derived by re- 
quiring the equivalence of the body- and global-frame 
metrics in the body's buffer region B A . More specifically, 
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one substitutes the expansions of the body-frame poten- 
tials (|3.8[) and the global- frame potentials (|3.17[) . along 
with the coordinate transformation (|3.13|) , into the tensor 
transformation law for the metric p,18[) . using (|3.4p to 
express the metrics in terms of the potentials in both co- 
ordinate systems. Matching coefficients of the resultant 
multipole expansions gives the moment transformation 
formulae. Having used the body-frame gauge conditions 
(|3.15p to eliminate the wordline-data functions cka and 
/3a, one finds that the transformation formulae involve 
only the various moments and the CoM worldlines z\. 
A more detailed account of the procedure is given in RF 
[HI Sec. V. 

In after presenting the various moment transformation 
formulae in Secs. lBllB2l and El we show in Sec. (IB 4p 
how to use them to arrive at the translational equations 
of motion for an iV-body system. 

Note that all the quantities appearing below (moments 
and worldlines) are treated as functions of the global time 
coordinate t (with the argument supressed). Any quan- 
tities originally defined as functions of the body-frame 



time coordinate sa (like the body-frame moments or the 
worldlines z A ) can be converted to functions of i by using 
the coordinate transformation (|3. 131) at y l A = 0, which 
gives 

s A = s° A (t) = t - c- 2 a A (sA)\ SA=t + 0(c- 4 ). (Bl) 

This change of variables does affect the forms of any equa- 
tions presented here (or elsewhere in the paper), as we 
use /(sa) and f(t) to denote the same quantity. 



1. Body- frame multipole moments — » global- frame 
multipole moments 

First, via the metric transformation law, the global- 
frame multipole moments M^ A and Z^ L A (defined by 
(|3.16|) 1 are given in terms of the body- frame moments 
M\ and S\ (defined by (03) ) and the CoM worldlines 

z A , by 



ryih 



M% + - 



2/ 2 + 5/ — 5 



2l 3 + 7l 2 + 161 + 7 j L 
(Z + l)(2/ + 3) UaMa 



2l 2 + 171 - 
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,J<ai 
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l + l 



Mf 



4(2Z - 1) 
2^ + 1 



rf A M{ <L - 1 8 ai>i 



4/ 



T+T A b 



+ 0{c-% 



e^< a 'S A 1>] + 0(c~ 2 ), 



(B2) 



(B3) 



Recall that v\ — z A and a A — z A ; the accelerations 
a A here (and anywhere they appear in lPN-order terms) 
may be replaced here with their Newtonian values from 
(|3.2ip . The global- frame monopole tidal moments G gj A 
appearing here (and needed only with Newtonian accu- 
racy here) can be expressed in terms of the body-frame 



multipole moments of bodies B ^ A and the bodies' 
worldlines by using Eq. (|2.14p . 

The specific (nonzero) instances of these formulae 
needed in the M1-M2-S2-Q2 system are as follows (with 
M l 2 3 ee QH and S\ = S 1 ): 



AL 



M t /3 2 M 2 
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2 U1 



n 3k Q 3k \ + 0(c- 4 ), 
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; v 2- — )+o(c ), 



r 5 
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Z g,2 = 

Z g,2 - 
yijk _ 



4Miv[ + 0(c~ 2 ), 
AM 2 v 2 + 0{c- 2 ), 

2Q 2 3 -2e llk S k + 0(c- 2 ) 
12 
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4viQi k - ^v a 2 Q a 2 <3 8 k>l + 0{c 



Here, r = \z 2 — Z\ \ and n % = (z 2 — z\)/r, as in (|3.23p . 
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2. Global-frame multipole moments — » 
global-frame tidal moments 

Next, one can express the global- frame tidal moments 
G^ A (t) and Y* L A (t) (def) a body A in terms of the global- 



frame multipole moments M£ A and Z* L A (def) of all the 
other bodies B ^ A and the all the bodies' worldlincs 
Zq. For this purpose, rather than working directly with 
Gg A , it is easier to work with the tensors F^ A defined 
by' 



g.ext 



1=0 
oo 



1=0 ' 
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2(2/ + 3)c : 



:J?Jx - z A )" L \ + 0(c~ 4 ) 



(B6) 



c.f. (|3.17[) . The tensor j£ A will not be needed in our 
calculations (and contains no extra information). Note 
that Fg A and G^ A agree at Newtonian order, 

F g L A = G L ^ A + 0(c- 2 ). 

The Newtonian part of G^ A is also given in Eq. (12. 14)) . 



By equating the Eq. (|B6|) and the external part of 



Eq. p,16p . one finds that F^ A is given by 
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where the tensors N^ A and P^ A are given by 



Ka = Ka 
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+ 0(c~ J ). 



Similarly, in the gravito-magnetic sector, one finds 
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For the M1-M2-S2-Q2 system, using the non-zero 
global-frame multipole moments from the last subsection, 
one must compute the I = 0, 1, 2, 3 (resp. Z = 0, 1) cases 
of dB7j) and the I = 0,1,2 (resp. I = 0) cases of (|B8j) 
for A = 2, B = 1 (resp. A = 1, B = 2). The deriva- 
tives appearing here can be easily expressed in terms of 



r = \z% — Z\\ and n l = (z\ ~ z \)/ r via the identities 
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<L> 



d L r 



r (2Z + 1)!! H+ 1 ' 
2T^T dL r + ^T (5(a!a! - ldL - 2) r- 



3. Global-frame tidal moments — > body-frame tidal 
moments 



Finally, the metric transformation law gives the body- 
frame tidal moments G A and H A (defined by l|3.8p ) in 
terms of the global-frame tidal moments F^ A and Y^ A 
(defined by (|B6|) . (|3.17|) V First, the gravito-magentic 
tidal moments can be found from 



ttL _ \rjk<L—l ai>jk 
n A — 1 A e > 



with the tensors Y A L given by 



Yf = (Sfk - - Av A «G« - UAf). (B9) 



Here, 8f K is the multi-Kronecker delta (5f K T^ K = T iL ), 
and 6f£ L> is the STF projector (S<^ L> T^ K = T <iL> ), 
and the nonzero 'inertial moments' A^ L are 
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Then, the gravito-electric tidal moments are given by 
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for I > 1 and by G A = for i = (c.f. ([SToc^ . The 
non-zero inertial moments needed here are 



AS 



A'-' 



a A + 
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G StA )a A + -v t ia A + 2G s , A v A 
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4. Translational equations of motion 



The results of Sees. IB 11 IB 21 and IB 31 allow one to ex- 
press the body-frame tidal moments G A and H A for a 
given body A in terms of the body-frame multipole mo- 
ments Mg and Sjf and CoM worldlines z l B of all bodies 



B in an A^-body system (eliminating all reference to the 
global frame moments). With this done, one can find the 
body's translational equation of motion, written only in 
terms of the Af J , S B , and z B , by using the law of motion 
(I3.20bp for the body- frame mass dipole M A . 



As the body-frame gauge condition (|3.15al) requires 
M A = (fixing the body's center of mass-energy to the 
body- frame origin), one procedes by setting the right- 
hand side of (I3.20b|) to zero. This yields an expres- 
sion for the acceleration a A = z\ of the body's 1PN- 
accurate global- frame CoM worldline z A (t). To see this 
more clearly, we can explicitly evaluate the I — case of 
the first term on the RHS of (|3.20b|) using the moment 
transformation formulae presented above. The result is 



M A a\ 



M A F„ 



1 1 

1=2 



Y 



S,A 



v A Y£ A + (2v\ - G StA )Gl A 



2 V A^g,A 



(v A + 3G g , A )c 



A A 



3G s , A v A 



(Bll) 

+ 0(c" 4 ), 



where g A represents the terms on the RHS of Eq. (13.20b)) 
except for the first. 

Appendix C: Derivation of system mass multipole 
moment formulae 

We derive here Eq. (|4.6[) . which gives the mass multi- 
pole moments of an A^-body system in terms of its 



global frame multipole moments M^ A and Z^ A . This 
parallels the derivation given in Sec. (IIV A[) for the mo- 
ments Zly S . We begin by using the Taylor series (|4.4[) to 
rewrite the global-frame scalar potential <£> g as given in 
Eq. (|4.1a|) in the form 



EE 

A l,k=0 
oo p 

EEE 

A p=0 fc=0 



(-1)'+* r 



l\k\ 

{-iy pi 

p\ kl{p-k)\ 



x\ 2c? 



d1{Ml A zld LK \x\) 



M <P-K K> 8 J_ , J_ pjl ( M (P-K K) 
M g,A Z A °P—\ + ^a°t \ M K,A Z A 



2c 2 



d P \x\ 



(CI) 



Here, To bring this into the form (|4.2a[) . which gives $ g in into its STF and trace parts. Using the identity 
terms of the system moments, we must decompose dp \ x | 1 w, „v 1 

d ijL \x\=d <ijL> \x\+ ^ + ^ + > S {l3 d L) - (C2) 
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and making the index change P — ¥ ijL, but only in the 
term resulting from the trace part of dp \x\, we find 



= - 



EEE 



A p=0k=0 

oo 1+2 



M 



<P-K„K> 



i (f + 2)i / ^fax-KK)\ ([+ 1)(/ + 2) i 

2c 2 2^ 2^ 2^ (/ + 2 )! fc!(Z + 2 - fc)! ' 1^ A ) 21 + 3 \x 

A 1=0 k=0 v 



<P>\X\ 



(C3) 



Though the sum over 2 should start at I = — 2 after _P — > 
r/L, the I = -1, -2 terms are killed by the (Z + 1)(Z + 2) 

factor. In the Ylk + =o t erm J the k indices K are to be 
chosen from the Z + 2 indices ij'L. This term can be 



simplified by explicitly performing the symmetrization 
over all Z+2 indices; using the fact that M^ A and <9l|:e| _1 
are STF, we have 



x 



2fc(Z + 2 - k) j<L -(K-l) K-l>j k(k-l) <L -(K-2)K-2>jj 

1V1„ A Z A + — — — TZM A Z A 

A (Z + 2)(Z + 1) S ' A A 



(Z + 2)(Z + 1) s- A 



x 



(C4) 



Using this identity, relabeling K — 1 — > K in the first term and K — 2 — > K in the second term, and adjusting 
summations appropriately, we find that the second line of (|C3I) can be written as 



1 



EEE 



(-iy 



2(2Z + 3)c 2 ^^^ Z! k\(l-k)\ 

y A 1=0 k=0 y ' 



(C5) 



Finally, we can compare (|C3[) (with the second line replaced by (|C5jl ) to the expression for $ g given in (|4.2a[) ; we see 
that the system's lPN-accurate mass multipoles must be given by 



Mi 



EE 



/! 



k\(l - k)\ 



A k=0 

1 2Z + 1 
_ ^(Z + l)(2Z + 3) 
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/i^ ys + 0( c - 4 ), 



as in Eq. (l4~6) . 



[1] C. Cutler, T. A. Apostolatos, L. Bildsten, L. S. Finn, 
E. E. Flanagan, D. Kennefick, D. M. Markovic, A. Ori, 
E. Poisson, G. J. Sussman, et al., Phys. Rev. Lett. 70, 
2984 (1993). 

[2] J. M. Lattimer and M. Prakash, Phys. Rep.-Rev. Sec. 

Phys. Lett. 442, 109 (2007). 
[3] M. D. Duez, Class. Quantum Gravity 27, 114002 (2010). 
[4] J. Faber, Class. Quantum Gravity 26, 114004 (2009). 
[5] J. A. Faber, P. Grandclement, F. A. Rasio, and 

K. Taniguchi, Phys. Rev. Lett. 89, 231102 (2002). 
[6] K. Kyutoku, M. Shibata, and K. Taniguchi, Phys. Rev. 

D 82, 044049 (2010). 



[7] M. Vallisneri, Phys. Rev. Lett. 84, 3519 (2000). 

[8] E. E. Flanagan and T. Hinderer, Phys. Rev. D 77, 021502 

(2008) . 

[9] T. Hinderer, B. D. Lackey, R. N. Lang, and J. S. Read, 
Phys. Rev. D 81, 123016 (2010). 
[10] T. Hinderer, Astrophys. J. 677, 1216 (2008). 
[11] T. Damour and A. Nagar, Phys. Rev. D 80, 084035 

(2009) . 

[12] T. Binnington and E. Poisson, Phys. Rev. D 80, 084018 
(2009). 

[13] T. Mora and C. M. Will, Phys. Rev. D 69, 104021 (2004). 
[14] T. Damour, M. Soffel, and C. Xu, Phys. Rev. D 43, 3273 



32 



(1991) . 

[15] T. Damour, M. Soffel, and C. Xu, Phys. Rev. D 45, 1017 

(1992) . 

[16] E. Racine and E. E. Flanagan, Phys. Rev. D 71, 044010 
(2005). 

[17] C. Xu, X. Wu, and G. Schafer, Phys. Rev. D 55, 528 
(1997). 

[18] X. J. Wu, Y. J. He, and C. M. Xu, Sci. China Ser. A- 

Math. Phys. Astron. 41, 1323 (1998). 
[19] J. S. Read, C. Markakis, M. Shibata, K. Uryu, J. D. E. 

Creighton, and J. L. Friedman, Phys. Rev. D 79, 124033 

(2009) . 

[20] L. Baiotti, T. Damour, B. Giacomazzo, A. Nagar, and 

L. Rezzolla, Phys. Rev. Lett. 105, 261101 (2010). 
[21] T. Damour and A. Nagar, Phys. Rev. D 81, 084016 

(2010) . 

[22] L. Blanchet, T. Damour, G. Esposito-Farese, and B. R. 



Iyer, Phys. Rev. Lett. 93, 091101 (2004). 
[23] L. Blanchet and T. Damour, Ann. Inst. Henri Poincarc- 

Phys. Thcor. 50, 377 (1989). 
[24] K. S. Thorne, Phys. Rev. D 58, 124031 (1998). 
[25] P. Purdue, Phys. Rev. D 60, 104054 (1999). 
[26] S. Weinberg, Gravitation and cosmology: Principles and 

applications of the general theory of relativity (1972). 
[27] K. S. Thorne and J. B. Hartle, Phys. Rev. D 31, 1815 

(1985). 

[28] V. Brumberg, Celest. Mech. Dyn. Astron. 99, 245 (2007). 
[29] T. Damour, C. R. Acad. Sci. Ser. Ii. 294, 1355 (1982). 
[30] S. Taylor and E. Poisson, Phys. Rev. D 78, 084016 
(2008). 

[31] A. Buonanno and T. Damour, Phys. Rev. D 59, 084006 
(1999). 



